Light field device, optical aberration compensation or simulation rendering method and vision testing system using same

ABSTRACT

Described are various embodiments of a light field device, optical aberration compensation or simulation rendering method and vision testing system using same. In one embodiment, the device comprises a digital display comprising a set of pixels; an array of light field shaping elements (LFSE) disposed relative to the set of pixels so to at least partially govern a light field emanated thereby; and a digital processor operable to: receive as input one or more higher order aberration parameters digitally defining a higher order aberration; for each given pixel, identify an adjusted image plane location corresponding thereto given a corresponding LFSE corresponding thereto and given said one or more higher order aberration parameters, and associate therewith an adjusted image value designated for the adjusted image location; and render for each said given pixel said adjusted image value associated therewith.

CROSS REFERENCE TO RELATED APPLICATIONS

This application is a divisional of U.S. patent application Ser. No. 17/510,297 filed Oct. 25, 2021, which is a continuation of PCT/IB2021/051868 filed Mar. 5, 2021, which is a continuation-in-part of International Application Serial No. PCT/IB2020/057887 filed Aug. 22, 2020 and of International Application No. PCT/US2020/058383 filed Oct. 30, 2020, each of which are a respective continuation-in-part of U.S. patent application Ser. No. 16/810,143 filed Mar. 5, 2020, now U.S. Pat. No. 10,761,604, which is a continuation-in-part of U.S. patent application Ser. No. 16/551,572 filed Aug. 26, 2019, now U.S. Pat. No. 10,636,116, a continuation-in-part of International Application No. PCT/IB2019/058955 filed Oct. 21, 2019, and a continuation-in-part of U.S. patent application Ser. No. 16/569,137 filed Sep. 12, 2019, now U.S. Pat. No. 10,642,355, the latter of which is a continuation of U.S. patent application Ser. No. 16/510,673 filed Jul. 12, 2019, now U.S. Pat. No. 10,474,235, which is a continuation of U.S. patent application Ser. No. 16/259,845 filed Jan. 28, 2019, now U.S. Pat. No. 10,394,322, and which claims priority to Canadian Patent Application No. 3,021,636 filed Oct. 22, 2018, the entire disclosure of each of which is hereby incorporated herein by reference.

U.S. patent application Ser. No. 17/510,297 is also a continuation-in-part of PCT/IB2020/057985 filed Aug. 26, 2020, which is a continuation of U.S. patent application Ser. No. 16/854,787 filed Apr. 21, 2020, now U.S. Pat. No. 10,860,099, the entire disclosure of each of which is also hereby incorporated herein by reference.

FIELD OF THE DISCLOSURE

The present disclosure relates to digital displays, and in particular, to a light field device, optical aberration compensation or simulation rendering method and vision testing system using same.

BACKGROUND

Individuals routinely wear corrective lenses to accommodate for reduced vision acuity in consuming images and/or information rendered, for example, on digital displays provided, for example, in day-to-day electronic devices such as smartphones, smart watches, electronic readers, tablets, laptop computers and the like, but also provided as part of vehicular dashboard displays and entertainment systems, to name a few examples. The use of bifocals or progressive corrective lenses is also commonplace for individuals suffering from near and farsightedness.

The operating systems of current electronic devices having graphical displays offer certain “Accessibility” features built into the software of the device to attempt to provide users with reduced vision the ability to read and view content on the electronic device. Specifically, current accessibility options include the ability to invert images, increase the image size, adjust brightness and contrast settings, bold text, view the device display only in grey, and for those with legal blindness, the use of speech technology. These techniques focus on the limited ability of software to manipulate display images through conventional image manipulation, with limited success.

The use of 4D light field displays with lenslet arrays or parallax barriers to correct visual aberrations have since been proposed by Pamplona et al. (PAMPLONA, V., OLIVEIRA, M., ALIAGA, D., AND RASKAR, R.2012. “Tailored displays to compensate for visual aberrations.” ACM Trans. Graph. (SIGGRAPH) 31.). Unfortunately, conventional light field displays as used by Pamplona et al. are subject to a spatio-angular resolution trade-off; that is, an increased angular resolution decreases the spatial resolution. Hence, the viewer sees a sharp image but at the expense of a significantly lower resolution than that of the screen. To mitigate this effect, Huang et al. (see, HUANG, F.-C., AND BARSKY, B. 2011. A framework for aberration compensated displays. Tech. Rep. UCB/EECS-2011-162, University of California, Berkeley, December; and HUANG, F.-C., LANMAN, D., BARSKY, B. A., AND RASKAR, R. 2012. Correcting for optical aberrations using multi layer displays. ACM Trans. Graph. (SiGGRAPH Asia) 31, 6, 185:1-185:12) proposed the use of multilayer display designs together with prefiltering. The combination of prefiltering and these particular optical setups, however, significantly reduces the contrast of the resulting image.

Finally, in U.S. Patent Application Publication No. 2016/0042501 and Fu-Chung Huang, Gordon Wetzstein, Brian A. Barsky, and Ramesh Raskar. “Eyeglasses-free Display: Towards Correcting Visual Aberrations with Computational Light Field Displays”. ACM Transaction on Graphics, 33(4), August 2014, the entire contents of each of which are hereby incorporated herein by reference, the combination of viewer-adaptive pre-filtering with off-the-shelf parallax barriers has been proposed to increase contrast and resolution, at the expense however, of computation time and power.

Optical devices, such as refractors and phoropters, are commonly used to test or evaluate the visual acuity of its users, for example, in the prescription of corrective eyewear, contact lenses or intraocular implants.

This background information is provided to reveal information believed by the applicant to be of possible relevance. No admission is necessarily intended, nor should be construed, that any of the preceding information constitutes prior art or forms part of the general common knowledge in the relevant art.

SUMMARY

The following presents a simplified summary of the general inventive concept(s) described herein to provide a basic understanding of some aspects of the disclosure. This summary is not an extensive overview of the disclosure. It is not intended to restrict key or critical elements of embodiments of the disclosure or to delineate their scope beyond that which is explicitly or implicitly described by the following description and claims.

A need exists for a light field device, optical aberration compensation or simulation rendering method and vision testing system using same that overcome some of the drawbacks of known techniques, or at least, provides a useful alternative thereto. For example, a need exists for solutions that can address or can be operated to test or simulate corrective measures for higher order aberrations or similar conditions, to name a few examples. Some aspects of disclosure provide embodiments of such systems, methods, and displays.

In accordance with one aspect, there is provided a computer-implemented method, automatically implemented by one or more digital processors, to automatically adjust user perception of an input image to be rendered via a set of pixels and a corresponding array of light field shaping elements (LFSE), wherein the user perception is affected by an optical aberration, the method comprising: for each given pixel of at least some of the set of pixels: digitally projecting an adjusted image ray trace from said given pixel towards a user pupil location, given a corresponding LFSE corresponding thereto, to intersect a virtual phase element digitally defined between said given pixel and the user's retina to virtually produce a designated optical phase shift corresponding to the optical aberration; digitally computing a corresponding refracted ray trace resulting from intersection of said adjusted image ray trace with said virtual phase element; digitally projecting said refracted ray trace to intersect an adjusted image plane location on a designated adjusted image plane; associating with said given pixel an adjusted image value designated for said adjusted image location; and rendering for each said given pixel said adjusted image value associated therewith.

In one embodiment, the adjusted image plane comprises a pixel plane.

In one embodiment, the optical aberration is defined by an eye focus parameter and a higher order aberration parameter, wherein said virtual phase element is defined as a function of said higher order aberration parameter, wherein said adjusted image plane comprises a virtual image plane virtually positioned relative to the digital display at a designated minimum viewing distance designated such that said perceptively adjusted version of the input image is adjusted to accommodate said eye focus parameter.

In one embodiment, the aberration is defined by an eye focus parameter and a higher order aberration parameter, wherein said virtual phase element is defined as a function of said higher order aberration parameter, wherein said adjusted image plane is designated as a user retinal plane, and wherein said refracted ray trace is projected to intersect with said user retinal plane by further redirecting said refracted image ray trace given said pupil location in accordance with said eye focus parameter.

In one embodiment, the virtual phase element is defined by a variable dioptric power, and wherein said refracted ray trace is refracted as a function of a given dioptric power defined for an intersection portion of said virtual phase element intersected by said adjusted image ray trace.

In one embodiment, the variable dioptric power is profiled in accordance with a higher order aberration of the user's eye.

In one embodiment, the variable dioptric power is profiled to simulate a progressive or bifocal lens.

In one embodiment, the virtual phase element is defined by an input surface and an output surface, and wherein said refracted ray trace is digitally refracted to output said by sequentially computing a first refraction event at said input surface and a second refraction event at said output surface.

In one embodiment, at least one of said input surface or said output surface has a variable profile.

In one embodiment, the virtual phase element is defined by a height function characterizing said variable profile.

In one embodiment, the height function is derived from a set of Zernike coefficients preset so as to reproduce said designated optical phase shift.

In one embodiment, the virtual phase element is further defined by a normal vector function derived from said height function that characterizes an orientation of an output surface of said virtual phase element with respect to any location on an input surface of said virtual phase element.

In one embodiment, the computing sequentially a first refraction event at said input surface and a second refraction event at said output surface comprises: identifying an input location of said adjusted ray trace on said input surface; computing a first refracted ray from said first refraction event at said input location; at least partially projecting said first refracted ray inside said virtual phase element; and computing said second refraction event using said normal vector function.

In one embodiment, the at least partially projecting said first refracted ray inside said virtual phase element comprises: tracing said first refracted ray until it reaches a location on said output surface; wherein said normal vector function in said computing said second refraction event is the normal vector function at the input surface location corresponding to the location on said output surface.

In one embodiment, the at least partially projecting said first refracted ray inside said phase element comprises: deriving a point on said output surface corresponding to said input location on said input surface but translated perpendicularly a value equal to said height function for that input location; wherein said normal vector function in said computing said second refraction event is the normal vector function at said input location.

In one embodiment, the at least partially projecting said first refracted ray inside said phase element comprises: partially tracing said first refracted ray a length equal to said height function for that input location to a partial location; and translating said first refracted ray perpendicularly to said input surface towards said output surface until it reaches said second surface; wherein said normal vector function in said computing said second refraction event is said normal vector function at the input surface location corresponding to said partial location on said output surface.

In one embodiment, the height function and said normal vector function are represented in memory in a data object comprising, respectively, a set of height values and associated normal vector values of said output surface corresponding to a plurality of surface locations on said input surface; wherein using said height function and said normal vector function involves fetching a height value and associated normal vector value corresponding to the surface location of said plurality of surface locations that is closest to said input location.

In one embodiment, the data object is a texture data object.

In one embodiment, the input image comprises distinct image portions, each image portion being associated with a distinct virtual phase element, and wherein said adjusted image ray trace for said distinct image portions intersects a corresponding one of said distinct virtual phase element to be processed accordingly.

In one embodiment, the virtual phase element comprises a plurality of virtual phase elements, and wherein said adjusted image ray trace is sequentially processed to intersect sequential ones of said virtual phase elements in its path.

In one embodiment, the optical aberration is defined by at least one of an astigmatism parameter, an eye focus parameter, or a higher order aberration parameter.

In one embodiment, the virtual phase element is virtually located at a location between said given pixel and a user pupil location so as to at least partially simulate an intervening corrective lens defined to produce said designated phase shift in at least partially addressing the optical aberration.

In one embodiment, the virtual phase element is virtually located so as to at least partially simulate an optical correction, at least partially defined by said designated phase shift when applied to compensate for the optical aberration, to be experienced by the user in selecting a corrective lens.

In one embodiment, the virtual phase element is virtually located at a user's eye pupil or lens location so to simulate an optical correction to be experienced by the user after undergoing a corresponding eye surgery or intraocular implant.

In accordance with another aspect, there is provided a computer-readable medium comprising instructions for implementation by one or more digital processors to automatically adjust user perception of an input image to be rendered via a set of pixels and a corresponding array of light field shaping elements (LFSE), wherein the user perception is affected by an optical aberration, by implementing a method a described above.

In accordance with another aspect, there is provided a device to automatically adjust user perception of an input image, wherein the user perception is affected by an optical aberration, the device comprising: a digital display comprising a set of pixels; an array of light field shaping elements (LFSE) disposed relative to said set of pixels so to at least partially govern a light field emanated thereby; a digital processor; and a computer-readable medium as defined above.

In accordance with another aspect, there is provided a digital display device to automatically adjust user perception of an input image, wherein the user perception is affected by an optical aberration, the device comprising: a digital display comprising a set of pixels; an array of light field shaping elements (LFSE) disposed relative to said set of pixels so to at least partially govern a light field emanated thereby; and a digital processor operable to: receive as input an optical aberration parameter; digitally define a virtual phase element based at least in part on said optical aberration parameter to be virtually disposed between said pixels and the user's retina to virtually produce a designated optical phase shift corresponding to the optical aberration; virtually trace a ray from each given pixel of at least some of said pixels to virtually interface with said virtual phase element so to define a corresponding refracted ray trace that intersects a respective adjusted image plane location; associate with said given pixel an adjusted image value designated for said respective adjusted image location; and render for each said given pixel said adjusted image value associated therewith.

In one embodiment, the optical parameter comprises one or more polynomial coefficients.

In one embodiment, the polynomial coefficients comprise Zernike coefficients.

In one embodiment, the virtual phase element is defined by a data object comprising a set of height values and associated normal vector values corresponding with a plurality of virtual phase element surface locations.

In one embodiment, the data object comprises texture data object.

In one embodiment, the optical aberration parameter comprises an eye focus parameter and a higher order aberration parameter, wherein said virtual phase element is defined, at least in part, as a function of both said eye focus parameter and said higher order aberration parameter.

In one embodiment, the optical aberration parameter comprises an eye focus parameter and a higher order aberration parameter, wherein said virtual phase element is defined, at least in part, as a function of said higher order aberration parameter, and wherein an adjusted image plane comprising each said adjusted image plane location is automatically designated as a function of said eye focus parameter so to distinctly compensate therefor.

In one embodiment, the virtual phase element is virtually located at a location between said given pixel and a user pupil location so as to at least partially simulate an intervening corrective lens defined to produce said designated phase shift in at least partially addressing the optical aberration.

In one embodiment, the virtual phase element is virtually located so as to at least partially simulate an optical correction, at least partially defined by said designated phase shift when applied to compensate for the optical aberration, to be experienced by the user in selecting a corrective lens.

In one embodiment, the virtual phase element is virtually located at a user's eye pupil or lens location so to simulate an optical correction to be experienced by the user after undergoing a corresponding eye surgery or intraocular implant.

In accordance with another aspect, there is provided a digital display device to automatically adjust user perception of an input image, wherein the user perception is affected by a higher order aberration, the device comprising: a digital display comprising a set of pixels; an array of light field shaping elements (LFSE) disposed relative to said set of pixels so to at least partially govern a light field emanated thereby; and a digital processor operable to: receive as input a set of higher order aberration parameters digitally defining the higher order aberration; for each given pixel of at least some of said pixels, identify an adjusted image plane location corresponding thereto given a corresponding LFSE corresponding thereto and given said set of higher order aberration parameters, and associate therewith an adjusted image value designated for said adjusted image location; and render for each said given pixel said adjusted image value associated therewith.

In one embodiment, the set of higher order aberration parameters comprises a set of Zernike coefficients.

In one embodiment, the user is further affected by a lower order aberration defined by an eye focus parameter, and wherein said adjusted image plane location is identified given said eye focus parameter and said set of higher order aberration parameters.

In one embodiment, the user is further affected by a lower order aberration defined by an eye focus parameter, and wherein an adjusted image plane comprising each said adjusted image plane location is automatically designated as a function of said eye focus parameter so to distinctly compensate therefor.

In one embodiment, the device simulates an intervening corrective lens that at least partially addresses the higher order aberration.

In one embodiment, the device simulates a higher order optical correction to be experienced by the user in selecting a corrective lens.

In one embodiment, the device simulates a higher order optical correction to be experience by the user after undergoing a corresponding eye surgery or intraocular implant.

Other aspects, features and/or advantages will become more apparent upon reading of the following non-restrictive description of specific embodiments thereof, given by way of example only with reference to the accompanying drawings.

BRIEF DESCRIPTION OF THE FIGURES

Several embodiments of the present disclosure will be provided, by way of examples only, with reference to the appended drawings, wherein:

FIGS. 1A and 1B are schematic diagrams of an exemplary light field vision testing or previewing system, in accordance with one embodiment;

FIGS. 2A to 2C schematically illustrate normal vision, blurred vision, and corrected vision in accordance with one embodiment, respectively;

FIGS. 3A and 3B are schematic diagrams of a light field display in which respective pixel subsets are aligned to emit light through a corresponding microlens or lenslet, in accordance with one embodiment;

FIGS. 4A to 4C are schematic diagrams of exemplary light field vision testing or previewing systems (e.g. refractors/phoropters), in accordance with different embodiments;

FIG. 5 is a plot of the angular resolution of an exemplary light field display as a function of the dioptric power generated, in accordance with one embodiment;

FIGS. 6A to 6D are schematic plots of the image quality generated by a light field refractor/phoropter as a function of the dioptric power generated by using in combination with the light field display (A) no refractive component, (B) one refractive component, (C) and (D) a multiplicity of refractive components;

FIGS. 7A, 7B and 7C are perspective views of exemplary light field refractors/phoropters, showing a casing thereof in cross-section (A and B) and a unit combining side-by-side two of the units (C) shown in 7A and 7B, in accordance with one embodiment;

FIG. 8 is a process flow diagram of an exemplary dynamic subjective vision testing method, in accordance with one embodiment;

FIG. 9 is a schematic diagram of an exemplary light field image showing two columns of optotypes at different dioptric power for the method of FIG. 8 , in accordance with one embodiment;

FIGS. 10A and 10B are process flow diagrams of exemplary input constant parameters and variables, respectively, for the ray-tracing rendering process of FIG. 11 , in accordance with one embodiment;

FIG. 11 is a process flow diagram of an illustrative ray-tracing rendering process, in accordance with one embodiment;

FIG. 12 is a process flow diagram illustrating a process step of FIG. 11 , in accordance with one embodiment;

FIG. 13 is a process flow diagram illustrating certain process steps of FIG. 11 , in accordance with one embodiment;

FIGS. 14A and 14B are schematic diagrams illustrating certain process steps of FIG. 11 , in accordance with one embodiment;

FIG. 15 is a schematic diagram illustrating the process steps of FIGS. 13 and 16 , in accordance with one embodiment;

FIG. 16 is a process flow diagram illustrating certain process steps of FIG. 11 , in accordance with one embodiment;

FIGS. 17A to 17D are schematic diagrams illustrating certain process steps of FIGS. 13 and 16 , in accordance with one embodiment;

FIGS. 18 and 19 are schematic diagrams illustrating ray-tracing in the context of non-parallel planes, in accordance with one embodiment;

FIGS. 20A and 20B are process flow diagrams of certain process steps of FIG. 11 for rendering a light field originating from multiple distinct virtual image planes, in accordance with one embodiment;

FIGS. 21A and 21B are process flow diagrams of certain process steps of FIGS. 20A and 20B for rendering a light field originating from multiple distinct virtual image planes, in accordance with one embodiment;

FIGS. 22A to 22D are schematic diagrams illustrating certain process steps of FIGS. 21A and 21B, in accordance with one embodiment;

FIGS. 23A and 23B are process flow diagrams of certain process steps of FIG. 11 for rendering a light field originating from multiple distinct focal planes, in accordance with one embodiment;

FIGS. 24A and 24B are schematic diagrams illustrating an example of a subjective visual acuity test using the ray-tracing rendering process of FIG. 23A or FIG. 23B, in accordance with one embodiment;

FIG. 25 is a schematic diagram of an exemplary virtual phase element for simulating optical aberrations, including higher order aberrations (HOAs), in accordance with one embodiment;

FIGS. 26A and 26B are schematic diagrams illustrating how a phase element of FIG. 25 may be used to simulate aberrations, including HOAs, in the eye (26A) or in a phase element located in front of the eye (20B), in accordance with one embodiment;

FIGS. 27A to 27D are schematic diagrams illustrating different exemplary shapes or profiles of the phase element of FIG. 25 , in accordance with one embodiment;

FIG. 28 is a schematic diagram illustrating a side-view of an exemplary phase element comprising a variable profile defined via a height or thickness function, in accordance with one embodiment;

FIG. 29 is a process flow diagram of an illustrative ray-tracing rendering process using the phase element of FIG. 25 and operable to compensate for HOAs, in accordance with one embodiment;

FIG. 30 is process flow diagram illustrating certain process steps of FIG. 29 , in accordance with one embodiment;

FIGS. 31A and 31B are process flow diagrams illustrating certain process steps of FIG. 30 , in accordance with one embodiment;

FIGS. 32A and 32B are plots of different profiles comparison of normally incident beam for raytracing through phase elements of FIG. 25 compared to diffraction with the corresponding aberration phase added to the wave front pre-diffraction. Normalized intensity profiles at the detector/image plane are shown, in accordance with one embodiment;

FIGS. 33A and 33B are exemplary plots of a height or thickness function for a phase element simulating astigmatism, in accordance with one embodiment;

FIG. 34 is a schematic diagram illustrating certain process steps of FIGS. 35 and 38 , in accordance with one embodiment;

FIG. 35 is a process flow diagram illustrating certain process steps of FIG. 29 , in accordance with one embodiment;

FIGS. 36A to 36C are schematic diagrams illustrating certain process steps of FIG. 35 , in accordance with one embodiment;

FIGS. 37A to 37C are schematic diagrams illustrating a phase element defined as a pre-computed set of height values and normal vectors, in accordance with one embodiment;

FIG. 38 is a process flow diagram illustrating certain process steps of FIG. 29 , in accordance with one embodiment;

FIG. 39 is a process flow diagram of an illustrative ray-tracing rendering process for rendering simultaneously multiple images or optotypes compensating for different sets of HOAs, in accordance with one embodiment;

FIGS. 40A and 40B are process flow diagrams illustrating certain process steps of FIG. 39 , in accordance with one embodiment;

FIG. 41 is a process flow diagram illustrating certain process steps of FIG. 40A, in accordance with one embodiment;

FIG. 42 is a process flow diagram illustrating certain process steps of FIG. 40B, in accordance with one embodiment;

FIG. 43 is a process flow diagram illustrating certain process steps of FIG. 40B, in accordance with one embodiment;

FIG. 44 is a schematic diagram illustrating the use of multiple phase elements simultaneously, in accordance with one embodiment;

FIG. 45 is a process flow diagram illustrating an exemplary unified light field rendering method incorporating phase elements, in accordance with one embodiment; and

FIG. 46 is a schematic diagram illustrating certain process steps of the method of FIG. 45 , in accordance with one embodiment.

Elements in the several figures are illustrated for simplicity and clarity and have not necessarily been drawn to scale. For example, the dimensions of some of the elements in the figures may be emphasized relative to other elements for facilitating understanding of the various presently disclosed embodiments. Also, common, but well-understood elements that are useful or necessary in commercially feasible embodiments are often not depicted in order to facilitate a less obstructed view of these various embodiments of the present disclosure.

DETAILED DESCRIPTION

Various implementations and aspects of the specification will be described with reference to details discussed below. The following description and drawings are illustrative of the specification and are not to be construed as limiting the specification. Numerous specific details are described to provide a thorough understanding of various implementations of the present specification. However, in certain instances, well-known or conventional details are not described in order to provide a concise discussion of implementations of the present specification.

Various apparatuses and processes will be described below to provide examples of implementations of the system disclosed herein. No implementation described below limits any claimed implementation and any claimed implementations may cover processes or apparatuses that differ from those described below. The claimed implementations are not limited to apparatuses or processes having all of the features of any one apparatus or process described below or to features common to multiple or all of the apparatuses or processes described below. It is possible that an apparatus or process described below is not an implementation of any claimed subject matter.

Furthermore, numerous specific details are set forth in order to provide a thorough understanding of the implementations described herein. However, it will be understood by those skilled in the relevant arts that the implementations described herein may be practiced without these specific details. In other instances, well-known methods, procedures and components have not been described in detail so as not to obscure the implementations described herein.

In this specification, elements may be described as “configured to” perform one or more functions or “configured for” such functions. In general, an element that is configured to perform or configured for performing a function is enabled to perform the function, or is suitable for performing the function, or is adapted to perform the function, or is operable to perform the function, or is otherwise capable of performing the function.

It is understood that for the purpose of this specification, language of “at least one of X, Y, and Z” and “one or more of X, Y and Z” may be construed as X only, Y only, Z only, or any combination of two or more items X, Y, and Z (e.g., XYZ, XY, YZ, ZZ, and the like). Similar logic may be applied for two or more items in any occurrence of “at least one . . . ” and “one or more . . . ” language.

Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.

Throughout the specification and claims, the following terms take the meanings explicitly associated herein, unless the context clearly dictates otherwise. The phrase “in one of the embodiments” or “in at least one of the various embodiments” as used herein does not necessarily refer to the same embodiment, though it may. Furthermore, the phrase “in another embodiment” or “in some embodiments” as used herein does not necessarily refer to a different embodiment, although it may. Thus, as described below, various embodiments may be readily combined, without departing from the scope or spirit of the innovations disclosed herein.

In addition, as used herein, the term “or” is an inclusive “or” operator, and is equivalent to the term “and/or,” unless the context clearly dictates otherwise. The term “based on” is not exclusive and allows for being based on additional factors not described, unless the context clearly dictates otherwise. In addition, throughout the specification, the meaning of “a,” “an,” and “the” include plural references. The meaning of “in” includes “in” and “on.”

As used in the specification and claims, the singular forms “a”, “an” and “the” include plural references unless the context clearly dictates otherwise.

The term “comprising” as used herein will be understood to mean that the list following is non-exhaustive and may or may not include any other additional suitable items, for example one or more further feature(s), component(s) and/or element(s) as appropriate.

The methods and systems described herein are directed, according to different embodiments, to light field rendering methods and systems capable of compensating for, or simulating, the presence of one or more optical aberrations, including higher order aberrations (HOAs), in, or disposed in relation to, a user's eye so to present a corrected, simulating or compensating image thereto. The methods and systems discussed below, in accordance with different embodiments, make use of a “virtual” optical element or phase element specifically tailored or defined so as to correctly reproduce or simulate the presence of these one or more aberrations.

For instance, a higher order aberration, as compared to lower zero, first or second-order aberrations, can often be defined by one or more polynomial functions that, alone or in combination, can be used to define a refractive response resulting from the higher order aberration. Examples of higher order aberrations may include, but are not limited to, trefoil, coma, and spherical aberrations that may result in poor night vision, glare, halos, blurring, starburst patterns or double vision. In these or other embodiments, lower order aberrations, such as those common to astigmatism (e.g. defined by cylindrical and axis parameters), myopia or hyperopia (e.g. defined by limited or constrained eye focus parameters), may be addressed, simulated or compensated for using the embodiments described herein, and that, independently or in combination. As such, addressed aberrations may be decomposed into lower and/or higher order aberrations as defined by respective or combined aberration parameters, whereby some may be otherwise addressed by alternate methods (e.g. within the simplified context of lower order aberrations), whereas other aberrations may be defined by complementary aberration parameters (e.g. higher order aberration parameters, etc.) and specifically addressed by the methods and approaches described herein. As such, in combination, all aberration parameters may be effectively addressed or simulated. These and other such considerations will become further apparent to the skilled artisan upon reading of the following illustrative examples.

As will be further described below, the herein described methods and systems may be equally applied to vision correction and/or for vision testing applications. For example, devices, displays and methods described herein may allow a user's perception of one or more input images (or input image portions), where each image or image portion is virtually located or perceived to be at a distinct image plane/depth location, to be adjusted or altered using the described light field display technology, again allowing for corrective assessment and/or previewing capabilities in respect of this viewer's reduced visual acuity or aberrations.

Some of the herein described embodiments provide for digital display devices, or devices encompassing such displays, for use by users having reduced visual acuity, whereby images ultimately rendered by such devices can be dynamically processed to accommodate the user's reduced visual acuity so that they may more comfortably consume rendered images without the use of corrective eyewear, contact lenses, or surgical intervention, as would otherwise be required. As noted above, embodiments are not to be limited as such as the notions and solutions described herein may also be applied to other technologies in which a user's perception of an input image to be displayed can be altered or adjusted via the light field display. Again, similar implementation of the herein described embodiments can allow for implementation of digitally adaptive vision tests or corrective or adaptive vision previews or simulations such that individuals with such reduced visual acuity can be exposed to distinct perceptively adjusted versions of an input image(s) to subjectively ascertain a potentially required or preferred vision correction.

Generally, digital displays as considered herein will comprise a set of image rendering pixels and a corresponding set of light field shaping elements that at least partially govern a light field emanated thereby to produce a perceptively adjusted version of the input image, notably distinct perceptively adjusted portions of an input image or input scene, which may include distinct portions of a same image, a same 2.5D/3D scene, or distinct images (portions) associated with different image depths, effects and/or locations and assembled into a combined visual input. For simplicity, the following will generally consider distinctly addressed portions or segments as distinct portions of an input image, whether that input image comprises a singular image having distinctly characterized portions, a digital assembly of distinctly characterized images, overlays, backgrounds, foregrounds or the like, or any other such digital image combinations.

In some examples, light field shaping elements may take the form of a light field shaping layer or like array of optical elements to be disposed relative to the display pixels in at least partially governing the emanated light field. As described in further detail below, such light field shaping layer elements may take the form of a microlens and/or pinhole array, or other like arrays of optical elements, or again take the form of an underlying light field shaping layer, such as an underlying array of optical gratings or like optical elements operable to produce a directional pixelated output.

Within the context of a light field shaping layer, as described in further detail below in accordance with some embodiments, the light field shaping layer can be disposed at a pre-set distance from the pixelated display so to controllably shape or influence a light field emanating therefrom. For instance, each light field shaping layer can be defined by an array of optical elements centered over a corresponding subset of the display's pixel array to optically influence a light field emanating therefrom and thereby govern a projection thereof from the display medium toward the user, for instance, providing some control over how each pixel or pixel group will be viewed by the viewer's eye(s). As will be further detailed below, arrayed optical elements may include, but are not limited to, lenslets, microlenses or other such diffractive optical elements that together form, for example, a lenslet array; pinholes or like apertures or windows that together form, for example, a parallax or like barrier; concentrically patterned barriers, e.g. cut outs and/or windows, such as a to define a Fresnel zone plate or optical sieve, for example, and that together form a diffractive optical barrier (as described, for example, in Applicant's co-pending U.S. application Ser. No. 15/910,908, the entire contents of which are hereby incorporated herein by reference); and/or a combination thereof, such as for example, a lenslet array whose respective lenses or lenslets are partially shadowed or barriered around a periphery thereof so to combine the refractive properties of the lenslet with some of the advantages provided by a pinhole barrier.

In operation, the display device will also generally invoke a hardware processor operable on image pixel (or subpixel) data for an image to be displayed to output corrected or adjusted image pixel data to be rendered as a function of a stored characteristic of the light field shaping elements and/or layer (e.g. layer distance from display screen, distance between optical elements (pitch), absolute relative location of each pixel or subpixel to a corresponding optical element, properties of the optical elements (size, diffractive and/or refractive properties, etc.), or other such properties, and a selected vision correction or adjustment parameter related to the user's reduced visual acuity or intended viewing experience. While light field display characteristics will generally remain static for a given implementation (i.e. a given shaping element and/or layer will be used and set for each device irrespective of the user), image processing can, in some embodiments, be dynamically adjusted as a function of the user's visual acuity or intended application so to actively adjust a distance of a virtual image plane, or perceived image on the user's retinal plane given a quantified user eye focus or like optical aberration(s), induced upon rendering the corrected/adjusted image pixel data via the static optical layer and/or elements, for example, or otherwise actively adjust image processing parameters as may be considered, for example, when implementing a viewer-adaptive pre-filtering algorithm or like approach (e.g. compressive light field optimization), so to at least in part govern an image perceived by the user's eye(s) given pixel or subpixel-specific light visible thereby through the layer.

With reference to FIGS. 1A and 1B, and in accordance with different embodiments, an exemplary subjective vision testing device/system (interchangeably referred to as a corrective vision previewing device/system), generally referred to using the numeral 100, will now be described. At the heart of this system is a light field vision testing device such as a light field refractor or phoropter device 102. Generally, light field refractor 102 is a device comprising, as mentioned above, a light field display 104 and which is operable to display or generate one or more images, including optotypes, to a user or patient having his/her vision acuity (e.g. refractive error) tested.

In some embodiments, as illustrated in FIG. 1B, light field display 104 comprises a light field shaping layer (LFSL) 108 overlaid or placed in front of a digital pixel display 110 (i.e. LCD, LED, OLED, etc.). For the sake of illustration, the following embodiments will be described within the context of a LFSL 108 defined, at least in part, by a lenslet array comprising an array of microlenses (also interchangeably referred to herein as lenslets) that are each disposed at a distance from a corresponding subset of image rendering pixels in an underlying digital display. It will be appreciated that while a light field shaping layer may be manufactured and disposed as a digital screen overlay, other integrated concepts may also be considered, for example, where light field shaping elements are integrally formed or manufactured within a digital screen's integral components such as a textured or masked glass plate, beam-shaping light sources (e.g. directional light sources and/or backlit integrated optical grating array) or like component.

Accordingly, each lenslet will predictively shape light emanating from these pixel subsets to at least partially govern light rays being projected toward the user by the display device. As noted above, other light field shaping layers may also be considered herein without departing from the general scope and nature of the present disclosure, whereby light field shaping will be understood by the person of ordinary skill in the art to reference measures by which light, that would otherwise emanate indiscriminately (i.e. isotropically) from each pixel group, is deliberately controlled to define predictable light rays that can be traced between the user and the device's pixels through the shaping layer.

For greater clarity, a light field is generally defined as a vector function that describes the amount of light flowing in every direction through every point in space. In other words, anything that produces or reflects light has an associated light field. The embodiments described herein produce light fields from an object that are not “natural” vector functions one would expect to observe from that object. This gives it the ability to emulate the “natural” light fields of objects that do not physically exist, such as a virtual display located far behind the light field display.

In one example, to apply this technology to vision correction, consider first the normal ability of the lens in an eye, as schematically illustrated in FIG. 2A, where, for normal vision, the image is to the right of the eye (C) and is projected through the lens (B) to the retina at the back of the eye (A). As comparatively shown in FIG. 2B, the poor lens shape and inability to accommodate (F) in presbyopia causes the image to be focused past the retina (D) forming a blurry image on the retina (E). The dotted lines outline the path of a beam of light (G). Naturally, other optical aberrations present in the eye will have different impacts on image formation on the retina. To address these aberrations, a light field display 104, in accordance with some embodiments, projects the correct sharp image (H) on the retina for an eye with a crystalline lens which otherwise could not accommodate sufficiently to produce a sharp image. The other two light field pixels (I) and (J) are drawn lightly, but would otherwise fill out the rest of the image.

As will be appreciated by the skilled artisan, a light field as seen in FIG. 2C cannot be produced with a ‘normal’ two-dimensional display because the pixels' light field emits light isotopically. Instead it is necessary to exercise tight control on the angle and origin of the light emitted, for example, using a microlens array or other light field shaping layer such as a parallax barrier, or combination thereof. Following with the example of a microlens array for LFSL 106, FIG. 3A schematically illustrates a single light field pixel defined by a convex microlens 302 disposed at its focus from a corresponding subset of pixels in a digital pixel display 108 to produce a substantially collimated beam of light emitted by these pixels, whereby the direction of the beam is controlled by the location of the pixel(s) relative to the microlens. The single light field pixel produces a beam similar to that shown in FIG. 2C where the outside rays are lighter and the majority inside rays are darker. The digital pixel display 108 emits light which hits the microlens 302 and it results in a beam of substantially collimated light (A).

Accordingly, upon predictably aligning a particular microlens array with a pixel array, a designated “circle” of pixels will correspond with each microlens and be responsible for delivering light to the pupil through that lens. FIG. 3B schematically illustrates an example of a light field display assembly in which a LFSL 106 sits above a pixel display 108 to have pixels 304 emit light through the microlens array. A ray-tracing algorithm can thus be used to produce a pattern to be displayed on the pixel array below the microlens in order to create the desired virtual image that will effectively correct for the viewer's reduced visual acuity.

As will be detailed further below, the separation between the LFSL 106 and the pixel array 108 as well as the pitch of the lenses can be selected as a function of various operating characteristics, such as the normal or average operating distance of the display, and/or normal or average operating ambient light levels.

In some embodiments, LFSL 106 may be a microlens array (MLA) defined by a hexagonal array of microlenses or lenslet disposed so to overlay a corresponding square pixel array of digital pixel display 108. In doing so, while each microlens can be aligned with a designated subset of pixels to produce light field pixels as described above, the hexagonal-to-square array mismatch can alleviate certain periodic optical artifacts that may otherwise be manifested given the periodic nature of the optical elements and principles being relied upon to produce the desired optical image corrections. Conversely, a square microlens array may be favoured when operating a digital display comprising a hexagonal pixel array.

In some embodiments, the MLA may further or alternatively be overlaid or disposed at an angle (rotation) relative to the underlying pixel array, which can further or alternatively alleviate period optical artifacts.

In yet some further or alternative embodiments, a pitch ratio between the microlens array and pixel array may be deliberately selected to further or alternatively alleviate periodic optical artifacts. For example, a perfectly matched pitch ratio (i.e. an exact integer number of display pixels per microlens) is most likely to induce periodic optical artifacts, whereas a pitch ratio mismatch can help reduce such occurrences.

Accordingly, in some embodiments, the pitch ratio will be selected to define an irrational number, or at least, an irregular ratio, so to minimize periodic optical artifacts. For instance, a structural periodicity can be defined so to reduce the number of periodic occurrences within the dimensions of the display screen at hand, e.g. ideally selected so to define a structural period that is greater than the size of the display screen being used.

While this example is provided within the context of a microlens array, similar structural design considerations may be applied within the context of a parallax barrier, diffractive barrier or combination thereof. In some embodiments, light field display 104 can render dynamic images at over 30 frames per second on the hardware in a smartphone.

Accordingly, a display device as described above and further exemplified below, can be configured to render a corrected or adjusted image via the light field shaping layer that accommodates, tests or simulates for the user's visual acuity. By adjusting the image correction in accordance with the user's actual predefined, set or selected visual acuity level, different users and visual acuity may be accommodated using a same device configuration, whereas adjusting such parameters for a given user may allow for testing for or simulation of different corrective or visual adjustment solutions. For example, by adjusting corrective image pixel data to dynamically adjust a virtual image distance below/above the display as rendered via the light field shaping layer, different visual acuity levels may be accommodated, and that, for an image input as a whole, for distinctly various portions thereof, or again progressively across a particular input.

As noted in the examples below, in some embodiments, light field rendering may be adjusted to effectively generate a virtual image on a virtual image plane that is set at a designated distance from an input user pupil location, for example, so to effectively push back, or move forward, a perceived image, or portion thereof, relative to the light field refractor device 102. In yet other embodiments, light field rendering may rather or alternatively seek to map the input image on a retinal plane of the user, taking into account visual aberrations, so to adaptively adjust rendering of the input image on the display device to produce the mapped effect. Namely, where the unadjusted input image would otherwise typically come into focus in front of or behind the retinal plane (and/or be subject to other optical aberrations), this approach allows to map the intended image on the retinal plane and work therefrom to address designated optical aberrations accordingly. Using this approach, the device may further computationally interpret and compute virtual image distances tending toward infinity, for example, for extreme cases of presbyopia. This approach may also more readily allow, as will be appreciated by the below description, for adaptability to other visual aberrations that may not be as readily modeled using a virtual image and image plane implementation. In both of these examples, and like embodiments, the input image is digitally mapped to an adjusted image plane (e.g. virtual image plane or retinal plane) designated to provide the user with a designated image perception adjustment that at least partially addresses designated visual aberrations. Naturally, while visual aberrations may be addressed using these approaches, other visual effects may also be implemented using similar techniques.

As an example of the effectiveness of the light field display in generating a diopter displacement (e.g. simulate the effect of looking through an optical component (i.e. a lens) of a given diopter strength or power) is shown in FIG. 5 , where a plot is shown of the angular resolution (in arcminutes) of an exemplary light field display comprising a 1500 ppi digital pixel display, as a function of the dioptric power of the light field image (in diopters). From this plot, it is clear that, in this particular example, the light field display is able to generate displacements (line 502) in diopters that have higher resolution corresponding to 20/20 vision (line 504) or better (e.g. 20/15—line 506) and close to (20/10—line 508)), here within a dioptric power range of 2 to 2.5 diopters.

Thus, in the context of a refractor 102, light field display 104 (in conjunction with the light field rendering or ray-tracing methods described below) may, according to different embodiments, be used to replace, at least in part, traditional optical components.

In some embodiments, the light field display can display a virtual image at optical infinity, meaning that any level of accommodation-based presbyopia (e.g. first order) can be corrected for. In some further embodiments, the light field display can both push the image back or forward, thus allowing for selective image corrections for both hyperopia (far-sightedness) and myopia (nearsightedness). In yet further embodiments as described below, variable displacements and/or accommodations may be applied as a function of non-uniform visual aberrations, or again to provide perceptive previewing or simulation of non-uniform or otherwise variable corrective powers/measures across a particular input or field of view.

However, the light field rendering system introduced above and the ray-tracing methods described below may also be used with other devices which may similarly comprise a light field display. For example, this may include a smartphone, tablets, e-readers, watches, televisions, GPS devices, laptops, desktop computer monitors, televisions, smart televisions, handheld video game consoles and controllers, vehicular dashboard and/or entertainment displays, and the like, without limitation.

Accordingly, any of the light field processing or ray-tracing methods described below, any modification thereto also discussed below, and related light field display solutions, can be equally applied to image perception adjustment solutions for visual media consumption, as they can for subjective vision testing solutions, or other technologically related fields of endeavour. As alluded to above, the light field display and rendering/ray-tracing methods discussed above may all be used to implement, according to various embodiments, a subjective vision testing device or system such as a phoropter or refractor. Indeed, a light field display may replace, at least in part, the various refractive optical components usually present in such a device. Thus, the vision correction light field ray tracing methods discussed below may equally be applied to render optotypes at different dioptric power or refractive correction by generating vision correction for hyperopia (far-sightedness) and myopia (nearsightedness), as was described above in the general case of a vision correction display. Light field systems and methods described herein, according to some embodiments, may be applied to create the same capabilities as a traditional instrument and to open a spectrum of new features, all while improving upon many other operating aspects of the device. For example, the digital nature of the light field display enables continuous changes in dioptric power compared to the discrete change caused by switching or changing a lens or similar; displaying two or more different dioptric corrections seamlessly at the same time; and, in some embodiments, the possibility of measuring higher-order aberrations and/or to simulate them for different purposes such as, deciding for free-form lenses, cataract surgery operation protocols, IOL choice, etc.

Going back to FIG. 1A, as producing a light field with angular resolution sufficient for accommodation correction over the full viewing ‘zone’ of a display would generally require an astronomically high pixel density, instead, a correct light field can be produced, in some embodiments, only at or around the location of the user's pupil(s). To do so, where a location or position of the user's eye is not otherwise rigidly constrained (e.g. within the context of a subject eye testing device or the like) the light field display can be paired with pupil tracking technology, as will be discussed below, to track a location of the user's eyes/pupils relative to the display. The display can then compensate for the user's eye location and produce the correct virtual image, for example, in real-time. Thus, in some embodiments, light field refractor 102 may include, integrated therein or interfacing therewith, a pupil/eye tracking system 110 to improve or enhance corrective image rendering by tracking a location of the user's eye(s)/pupil(s) (e.g. both or one, e.g. dominant, eye(s)) and adjusting light field corrections accordingly one or more eye/pupil tracking light sources, such as one or more infrared (IR) or near-IR (NIR) light source(s) to accommodate operation in limited ambient light conditions, leverage retinal retro-reflections, invoke corneal reflection, and/or other such considerations. For instance, different IR/NIR pupil tracking techniques may employ one or more (e.g. arrayed) directed or broad illumination light sources to stimulate retinal retro-reflection and/or corneal reflection in identifying a tracking a pupil location. Other techniques may employ ambient or IR/NIR light-based machine vision and facial recognition techniques to otherwise locate and track the user's eye(s)/pupil(s). To do so, one or more corresponding (e.g. visible, IR/NIR) cameras may be deployed to capture eye/pupil tracking signals that can be processed, using various image/sensor data processing techniques, to map a 3D location of the user's eye(s)/pupil(s). As mentioned above, in some embodiments, such eye/pupil tracking hardware/software may be integral to device 102, for instance, operating in concert with integrated components such as one or more front facing camera(s), onboard IR/NIR light source(s) (not shown) and the like. In other user environments, such as in a vehicular environment, eye/pupil tracking hardware may be further distributed within the environment, such as dash, console, ceiling, windshield, mirror or similarly-mounted camera(s), light sources, etc.

In one embodiment and as illustrated in FIG. 4A, light field refractor 102 may be configured with light field display 104 located relatively far away (e.g. one or more meters) from the user's eye currently being diagnosed. Note that the pointed line in FIGS. 4A to 4C is used to schematically illustrate the direction of the light rays emitted by light field display 104. Also illustrated is eye-tracker 110, which may be provided as a physically separate element, for example, installed in at a given location in a room or similar. In some embodiments, the noted eye/pupil tracker 110 may include the projection of IR markers/patterns to help align the patient's eye with the light field display. In some embodiments, a tolerance window (e.g. “eye box”) may be considered to limit the need to refresh the ray-tracing iteration. An exemplary value of the size of the eye box, in some embodiments, is around 6 mm, though smaller (e.g. 4 mm) or larger eye boxes may alternatively be set to impact image quality, stability or like operational parameters.

Going back to FIG. 1A, light field refractor 102 may also comprise, according to different embodiments and as will be further discussed below, one or more refractive optical components 112, a processing unit 114, a data storage unit or internal memory 116, one or more cameras 118, a power source 120, a network interface 122 for communicating via network to a remote database or server 124.

In some embodiments, power source 120 may comprise, for example, a rechargeable Li-ion battery or similar. In some embodiments, it may comprise an additional external power source, such as, for example, a USB-C external power supply. It may also comprise a visual indicator (screen or display) for communicating the device's power status, for example whether the device is on/off or recharging.

In some embodiments, internal memory 116 may be any form of electronic storage, including a disk drive, optical drive, read-only memory, random-access memory, or flash memory, to name a few examples. In some embodiments, a library of chart patterns (Snellen charts, prescribed optotypes, forms, patterns, or other) may be located in internal memory 116 and/or retrievable from remote server 124 via network interface 122.

In some embodiments, one or more optical components 112 may be used in combination with the light field display 104, for example to shorten the size of refractor 102 and still offer an acceptable range in dioptric power. The general principle is schematically illustrated in the plots of FIGS. 6A to 6D. In these schematic plots, the image quality (e.g. inverse of the angular resolution, higher is better) at which optotypes are small enough to be useful for vision testing in this plot is above horizontal line 602 which represents typical 20/20 vision. FIG. 6A shows the plot for the light field display only, where we see the characteristic two peaks corresponding to the smallest resolvable point, one of which was plotted in FIG. 5 (here inverted and shown as a peak instead of a basin), and where each region above the line may cover a few diopters of dioptric power, according to some embodiments. While the dioptric range may, in some embodiments, be more limited than needed when relying only on the light field display, it is possible to shift this interval by adding one or more refractive optical components. This is shown in FIG. 6B where the regions above the line 602 is shifted to the left (negative diopters) by adding a single lens in the optical path.

Thus, by using a multiplicity of refractive optical components 112 or by alternating sequentially between different refractive components 112 of increasing or decreasing dioptric power, it is possible to shift the center of the light field diopter range to any required value, as shown in FIG. 6C, and thus the image quality may be kept above line 602 for any required dioptric power as shown in FIG. 6D. In some embodiments, a range of 30 diopters from +10 to −20 may be covered for example. In the case of one or more reels of lenses, the lens may be switched for a given larger dioptric power increment, and the light field display would be used to provide a finer continuous change to accurately pin-point the required total dioptric power required to compensate for the patient's reduced visual acuity. This would still result in light field refractor 102 having a reduced number of refractive optical components compared to the number of components needed in a traditional refractor, while drastically enhancing the overall fine-tuning ability of the device.

One example, according to one embodiment, of such a light field refractor 102 is schematically illustrated in FIG. 4B, wherein the light field display 104 (herein shown comprising LFSL 106 and digital pixel display 108) is combined with a multiplicity of refractive components 112 (herein illustrated as a reel of lenses as an example only). By changing the refractive component used in combination with light field display 104, a larger dioptric range may be covered. This may also provide means to reduce the dimension of device 102 as mentioned above, making it more portable, so that all its internal components may be encompassed into a shell, housing or casing 402. In some embodiments, light field refractor 102 may thus comprise a durable ABS housing that may be shock and harsh-environment resistant. In some embodiments, light field refractor 102 may also comprise a telescopic feel for fixed or portable usage; optional mounting brackets, and/or a carrying case (not shown). In some embodiments, all components may be internally protected and sealed from the elements.

In some embodiments, casing 402 may further comprise a head-rest or similar (not shown) to keep the user's head still and substantially in the same location, thus, in such examples, foregoing the general utility of a pupil tracker or similar techniques by substantially fixing a pupil location relative to this headrest.

In some embodiments, it may also be possible to further reduce the size of device 102 by adding, for example, a mirror or any device which may increase the optical path. This is illustrated in FIG. 4C where the length of the device was reduced by adding a mirror 404. This is shown schematically by the pointed arrow which illustrates the light being emitted from pixel display 108 travelling through LFSL 106 before being reflected back by mirror through refractive components 112 and ultimately hitting the eye.

The skilled technician will understand that different examples of refractive components 112 may be include, without limitation, one or more lenses, sometimes arranged in order of increasing dioptric power in one or more reels of lenses similar to what is typically found in traditional refractors/phoropters; an electrically controlled fluid lens; active Fresnel lens; and/or Spatial Light Modulators (SLM). In some embodiments, additional motors and/or actuators (not shown) may be used to operate refractive components 112. The motors/actuators may be communicatively linked to processing unit 114 and power source 120, and operate seamlessly with light field display 102 to provide the required dioptric power.

For example, FIGS. 7A and 7B show a perspective view of an exemplary light field phoropter 102 similar to the one schematically shown in FIG. 3B, but wherein the refractive component 112 is an electrically tunable liquid lens. Thus, in this particular embodiment, no mechanical or moving component are used, which may result in the device being more robust. In some embodiments, the electrically tunable lens may have a range of ±13 diopters.

In one illustrative embodiment, a 1000 dpi display is used with a MLA having a 65 mm focal distance and 1000 μm pitch with the user's eye located at a distance of about 26 cm. A similar embodiment uses the same MLA and user distance with a 3000 dpi display.

Other displays having resolutions including 750 dpi, 1000 dpi, 1500 dpi and 3000 dpi were also tested or used, as were MLAs with a focal distance and pitch of 65 mm and 1000 μm, 43 mm and 525 μm, 65 mm and 590 μm, 60 mm and 425 μm, 30 mm and 220 μm, and 60 mm and 425 μm, respectively, and user distances of 26 cm, 45 cm or 65 cm.

Going back to FIG. 1A, in some embodiments, eye-tracker 110 may further comprise a digital camera, in which case it may be used to further acquire images of the user's eye to provide further diagnostics, such as pupillary reflexes and responses during testing for example. In other embodiments, one or more additional cameras 118 may be used to acquire these images instead. In some embodiments, light field refractor 102 may comprise built-in stereoscopic tracking cameras.

In some embodiments, feedback and/or control of the vision test being administered by system 100 may be given via a control interface 126. In some embodiments, the control interface 126 may comprise a dedicated handheld controller-like device 128. This controller 128 may be connected via a cable or wirelessly, and may be used by the patient directly and/or by an operator like an eye professional. In some embodiments, both the patient and operator may have their own dedicated controller 128. In some embodiments, the controller may comprise digital buttons, analog thumbstick, dials, touch screens, and/or triggers.

In some embodiments, control interface 126 may comprise a digital screen or touch screen, either on refractor 102 itself or part of an external module (not shown). In other embodiments, control interface 126 may let on or more external remote devices (i.e. computer, laptop, tablet, smartphone, remote, etc.) control light field refractor 102 via network interface 122. For example, remote digital device 130 may be connected to light field refractor 102 via a cable (e.g. USB cable, etc.) or wirelessly (e.g. via Wi-Fi, Bluetooth or similar) and interface with light field refractor 102 via a dedicated application, software or website (not shown). Such a dedicated application may comprise a graphical user interface (GUI), and may also be communicatively linked to remote database 124.

In some embodiments, the user or patient may give feedback verbally and the operator may control the vision test as a function of that verbal feedback. In some embodiments, refractor 102 may comprise a microphone (not shown) to record the patient's verbal communications, either to communicate them to a remote operator via network interface 122 or to directly interact with the device (e.g. via speech recognition or similar).

Going back to FIG. 1A, processing unit 114 may be communicatively connected to data storage 116, eye tracker 110, light field display 104 and refractive components 112. Processing unit 114 may be responsible for rendering one or more images or optotypes via light field display 104 and, in some embodiments, jointly control refractive components 112 to achieve a required total change in dioptric power. It may also be operable to send and receive data to internal memory 116 or to/from remote database 124 via network interface 122.

In some embodiments, diagnostic data may be automatically transmitted/communicated to remote database 124 or remote digital device 130 via network interface 122 through the use of a wired or wireless network connection. The skilled artisan will understand that different means of connecting electronic devices may be considered herein, such as, but not limited to, Wi-Fi, Bluetooth, NFC, Cellular, 2G, 3G, 4G, 5G or similar. In some embodiments, the connection may be made via a connector cable (e.g. USB including microUSB, USB-C, Lightning connector, etc.). In some embodiments, remote digital device 130 may be located in a different room, building or city.

In some embodiments, two light field refractors 102 may be combined side-by-side to independently measure the visual acuity of both left and right eye at the same time. An example is shown in FIG. 7C, where two units corresponding to the embodiment of FIG. 7A or 7B (used as an example only) are placed side-by-side or fused into a single device.

In some embodiments, a dedicated application, software or website may provide integration with third party patient data software. In some embodiments, software required to operate and installed on refractor 102 may be updated on-the-fly via a network connection and/or be integrated with the patient's smartphone app for updates and reminders.

In some embodiments, the dedicated application, software or website may further provide a remote, real-time collaboration platform between an eye professional and user/patient, and/or between different eye professionals. This may include interaction between different participants via video chat, audio chat, text messages, etc.

In some embodiments, light field refractor 102 may be self-operated or operated by an optometrist, ophthalmologist or other certified eye-care professional. For example, in some embodiments, a user/patient may use refractor 102 in the comfort of his/her own home, in a store or a remote location.

With reference to FIG. 8 and in accordance with one exemplary embodiment, a dynamic subjective vision testing method using vision testing system 100, generally referred to using the numeral 800, will now be described. As mentioned above, the use of a light field display enables refractor 102 to provide more dynamic and/or more modular vision tests than what is generally possible with traditional refractors/phoropters. Generally, method 800 seeks to diagnose a patient's reduced visual acuity and produce therefrom, in some embodiments, an eye prescription or similar.

In some embodiments, eye prescription information may include, for each eye, one or more of: distant spherical, cylindrical and/or axis values, and/or a near (spherical) addition value.

In some embodiments, the eye prescription information may also include the date of the eye exam and the name of the eye professional that performed the eye exam. In some embodiments, the eye prescription information may also comprise a set of vision correction parameter(s), as will be further discussed below, for operating any vision correction light field displays using the systems and methods described below. In some embodiments, the eye prescription may be tied to a patient profile or similar, which may contain additional patient information such as a name, address or similar. The patient profile may also contain additional medical information about the user. All information or data (i.e. set of vision correction parameter(s), user profile data, etc.) may be kept on external database 124. Similarly, in some embodiments, the user's current vision correction parameter(s) may be actively stored and accessed from external database 124 operated within the context of a server-based vision correction subscription system or the like, and/or unlocked for local access via the client application post user authentication with the server-based system.

Refractor 102 being, in some embodiments, portable, a large range of environments may be chosen to deliver the vision test (home, eye practitioner's office, etc.). At the start, the patient's eye may be placed at the required location. This may be done by placing his/her head on a headrest or by placing the objective (i.e. eyepiece) on the eye to be diagnosed. As mentioned above, the vision test may be self-administered or partially self-administered by the patient. For example, the operator (e.g. eye professional or other) may have control over the type of test being delivered, and/or be the person who generates or helps generate therefrom an eye prescription, while the patient may enter inputs dynamically during the test (e.g. by choosing or selecting an optotype, etc.).

As will be discussed below, light field rendering methods described herein generally requires an accurate location of the patient's pupil center. Thus, at step 802, such a location is acquired. In some embodiments, such a pupil location may be acquired via eye tracker 110, either once, at intervals, or continuously. In other embodiments, the location may be derived from the device or system's dimension. For example, in some embodiments, the use a head-rest and/or an eye-piece or similar provides an indirect means of deriving the pupil location. In some embodiments, refractor 102 may be self-calibrating and not require any additional external configuration or manipulation from the patient or the practitioner before being operable to start a vision test.

At step 804, one or more optotypes is/are displayed to the patient, at one or more dioptric power (e.g. in sequence, side-by-side, or in a grid pattern/layout). The use of light field display 104 offers multiple possibilities regarding how the images/optotypes are presented, and at which dioptric power each may be rendered. The optotypes may be presented sequentially at different dioptric power, via one or more dioptric power increments. In some embodiments, the patient and/or operator may control the speed and size of the dioptric power increments.

In some embodiments, optotypes may also be presented, at least in part, simultaneously on the same image but rendered at a different dioptric power. For example, FIG. 9 shows an example of how different optotypes may be displayed to the patient but rendered with different dioptric powers simultaneously. These may be arranged in columns or in a table or similar. In FIG. 9 , we see two columns of three optotypes (K, S, V), varying in size, as they are perceived by a patient, each column being rendered at different degrees of refractive correction (e.g. dioptric power). In this specific example, the optotypes on the right are being perceived as blurrier than the optotypes on the left.

Thus, at step 806, the patient would communicate/verbalize this information to the operator or input/select via, for example, control interface 126 the left column as the one being clearer. Thus, in some embodiments, method 800 may be configured to implement dynamic testing functions that dynamically adjust one or more displayed optotype's dioptric power in real-time in response to a designated input, herein shown by the arrow going back from step 808 to step 804 in the case where at step 808, the user or patient communicates that the perceived optotypes are still blurry or similar. In the case of sequentially presented optotypes, the patient may indicate when the optotypes shown are clearer. In some embodiments, the patient may control the sequence of optotypes shown (going back and forth as needed in dioptric power), and the speed and increment at which these are presented, until he/she identifies the clearest optotype. In some embodiments, the patient may indicate which optotype or which group of optotypes is the clearest by moving an indicator icon or similar within the displayed image.

In some embodiments, the optotypes may be presented via a video feed or similar.

In some embodiments, when using a reel of lenses or similar (for refractive components 112), discontinuous changes in dioptric power may be unavoidable. For example, the reel of lenses may be used to provide a larger increment in dioptric power, as discussed above. Thus, step 804 may in this case comprise first displaying larger increments of dioptric power by changing lens as needed, and when the clearest or less blurry optotypes are identified, fine-tuning with continuous or smaller increments in dioptric power using the light field display. In the case of optotypes presented simultaneously, the refractive components 112 may act on all optotypes at the same time, and the change in dioptric power between them may be controlled only by the light display 104. In some embodiments, for example when using an electrically tunable fluid lens or similar, the change in dioptric power may be continuous.

In some embodiments, eye images may be recorded during steps 802 to 806 and analyzed to provide further diagnostics. For example, eye images may be compared to a bank or database of proprietary eye exam images and analyzed, for example via an artificial intelligence (AI) or Machine-learning (ML) system or similar. This analysis may be done by refractor 102 locally or via a remote server or database 124.

Once the correct dioptric power needed to correct for the patient's reduced visual acuity is defined at step 810, an eye prescription or vision correction parameter(s) may be derived from the total dioptric power used to display the best perceived optotypes.

In some embodiments, the patient, an optometrist or other eye-care professional may be able to transfer the patient's eye prescription directly and securely to his/her user profile store on said server or database 124. This may be done via a secure website, for example, so that the new prescription information is automatically uploaded to the secure user profile on remote database 124. In some embodiments, the eye prescription may be sent remotely to a lens specialist or similar to have prescription glasses prepared.

In some embodiments, vision testing system 100 may also or alternatively be used to simulate compensation for other or higher-order aberrations. Indeed, the light field rendering methods described above may be used to compensation for higher order aberrations (HOA), and thus be used to validate externally measured or tested HOA via method 3600, in that a measured, estimated or predicted HOA can be dynamically compensated for using the system described herein and thus subjectively visually validated by the viewer in confirming whether the applied HOA correction satisfactorily addresses otherwise experienced vision deficiencies.

With reference to FIGS. 10A to 17D, and in accordance with different embodiments, an exemplary computationally implemented ray-tracing method for rendering an adjusted image via an array of light field shaping elements, in this example provided by a light field shaping layer (LFSL) disposed relative to a set of underlying display pixels, that accommodates or compensates for the user's reduced visual acuity will now be described. In this first example, for illustrative purposes, adjustment of a single image (i.e. the image as whole) is being implemented without consideration for distinct image portions and taking only into account spherical dioptric power changes. Further examples below will specifically address modification of the following example for adaptively adjusting distinct image portions and/or adjustments considering other aberrations such as HOA, and/or astigmatism.

In this exemplary embodiment and as shown in FIG. 10A, a set of constant parameters 1002 used for the light field rendering process may be pre-determined. These may include, for example, any data or parameters that are not expected to significantly change during a user's viewing session, between different viewing sessions or even between users, for instance. These may generally be based on the physical and functional characteristics of light field display 104 for which the method is to be implemented, as will be explained below. Similarly, as shown in FIG. 10B, every iteration of the rendering algorithm (i.e. when rendering a full light field image frame) may also use a set of input variables 1004 which are expected to change either at each rendering iteration (i.e. between frames) or at least between each user's viewing session.

As illustrated in FIG. 10A, the list of constant parameters 1002 may include, without limitations, the distance 1006 between pixel display 108 and the LFSL 106, the in-plane rotation angle 1008 between the frames of reference digital pixel display 108 and LFSL 106, the resolution 1010 and/or size 1011 of digital pixel display 108, the size 1012 of each individual pixel, the optical geometry 1014 of LFSL 106, the size or pitch of individual optical elements or units 1016 within LFSL 106 and optionally the subpixel layout 1018 of pixel display 108. Moreover, both the resolution 1010 and the size of each individual pixel 1012 of pixel display 108 may be used to pre-determine both the absolute size of the display in real units (i.e. in mm) and the three-dimensional position/location of each pixel within the display. In some embodiments where the subpixel layout 1018 is available, the position/location within pixel display 108 of each subpixel may also be pre-determined. These three-dimensional location/positions are usually calculated using a given frame of reference located somewhere within the plane of pixel display 108, for example a corner or the middle of the display, although other reference points may be chosen. Concerning the optical layer geometry 1014, different geometries may be considered, for example a hexagonal geometry. Finally, by combining the distance 1006, the rotation angle 1008, and the geometry 1014 with the optical element size 1016, it is possible to similarly pre-determine the three-dimensional location/position of each optical element of LFSL 106 with respect to the frame of reference of pixel display 108.

FIG. 10B meanwhile illustratively lists an exemplary set of input variables 1004 that may be used for the light field rendering methods described below, and may include any input data that may reasonably change during a user's single viewing session, for example between each image frame, or between different viewing sessions or even between different users (e.g. not related to the hardware specifications of light field display 104 or refractor 102). These may be automatically acquired by device 102 during normal operation or be inputted (for example by a user or person operating the device for the user) into the device as required. They may thus include without limitation: the image(s)/optotype(s) to be displayed 1020 (e.g. comprising pixel data such as on/off, colour, brightness, etc.), a three-dimensional pupil center location 1022 (e.g. in embodiments implementing active eye/pupil tracking methods) and/or a pupil size 1024.

The methods described below, according to different embodiments, also utilize vision correction parameters in the form of dioptric power so as to modulate the strength and nature of the compensation/correction generated by the light field image. These may include a spherical dioptric power 1026 (which may be derived indirectly, for example, from a minimum reading distance value 1028 as will be discussed below), but also, in some embodiments, one or more sets of a cylindrical power 1030 and a corresponding cylindrical axis angle 1032. In some embodiments, these input variables (spherical dioptric power 1026, cylindrical dioptric power 1026 and cylinder axis angle 1030) mirror the SPHERE, CYL and AXIS parameters used in a typical eye examination.

Moreover, in some embodiments, an eye depth value 1034 may also be used, either as an average value or customized for an individual user. In some embodiments, input image 1020, may be representative of one or more digital images to be displayed with digital pixel display 108. In addition, input image 1020 may generally be encoded in any data format used to store digital images known in the art.

Pupil center location 1022, in one embodiment, is the three-dimensional coordinates of at least one the user's pupils' center with respect to a given reference frame, for example a point on device 102 or digital pixel display 108 and may be derived from any eye/pupil tracking method known in the art via eye/pupil tracker 110. In some embodiments, the pupil center location 1022 may be determined prior to any new iteration of the rendering algorithm, or in other cases, at a lower framerate (and thus re-use the same location/position for two or more subsequent image frames). In some embodiments, only the pupil center location of a single user's eye may be determined, for example the user's dominant eye (i.e. the one that is primarily relied upon by the user). In some embodiments, this location/position, and particularly the associated pupil distance to the screen may otherwise or additionally be rather approximated or adjusted based on other contextual or environmental parameters, such as an average or preset user distance to the screen (e.g. typical reading distance for a given user or group of users; stored, set or adjustable driver distance in a vehicular environment; etc.).

In the illustrated embodiment, the minimum reading distance 1028 is defined as the minimal focus distance for reading that the user's eye(s) may be able to accommodate (i.e. able to view without discomfort). In some embodiments, different values of the minimum reading distance 1028 associated with different users may be entered, for example, as can other adaptive vision correction parameters be considered depending on the application at hand and vision correction being addressed. As mentioned above, in some embodiments, minimum reading distance 1028 may be derived from an eye prescription (e.g. glasses prescription or contact prescription) or similar. It may, for example, correspond to the near point distance corresponding to the uncorrected user's eye, which can be calculated from the prescribed corrective lens power assuming that the targeted near point was at 25 cm.

In some embodiments, as will be further described below, input variables 1004 may further comprise a set of Zernike coefficients 1036. These are usually used to define or characterize a set of HOAs to be corrected, as will be described below, in accordance with different embodiments. These may be chosen so as to describe one or more optical aberrations via a Zernike polynomial. In some embodiments, these may be used in addition to or in replacement of previously discussed variables such as spherical power 1026, cylindrical power 1030 or cylinder axis angle 1032. In other embodiments, at least some of the Zernike coefficients 1036 may be derived from spherical power 1026, cylindrical power 1030 or cylinder axis angle 1032. For example, in some embodiments, Zernike coefficients 1036 may be used to model HOA only, while spherical power may be simulated via the method 1100 described below. In some embodiments, in the exemplary case of astigmatism, the corresponding Zernike coefficients 1036 may be computed from spherical power 1026, cylindrical power 1030 or cylinder axis angle 1032, as will be shown below.

In addition, the methods for treating HOAs discussed below, in accordance with different embodiments, will make use of an optical or phase element to simulate the effects of HOAs. This element may have additional variables associated with it, for example at least one refractive index 1038 and a location/orientation 1040. More details regarding the use of this phase element (or a plurality of elements) will be given further below.

With added reference to FIGS. 11 to 16D, parameters 1002 and variables 1004 may be used in the light field ray-tracing method 1100, herein presented in accordance with different embodiments. While method 1100 and the steps described in FIG. 12 apply equally to both embodiments, steps 1102 and 1112 may be applied using either virtual image planes “virtually” located behind the display or on the retinal plane (via an eye focal plane (i.e. inside the eye)).

Moreover, method 1100 as illustrated in FIGS. 11 to 17D is designed to correct for spherical aberrations. However, other methods expanding on method 1100, for example that may also correct for astigmatism, will be described further below.

Moreover, for illustrative purposes, in this example, adjustment of a single image (i.e. the image as whole) is being implemented without consideration for distinct image portions. Further examples below will specifically address modification of the following example for adaptively adjusting distinct image portions.

In the illustrated embodiment of FIG. 11 , method 1100 begins with step 1102, in which the image to be displayed is pre-processed for subsequent ray-tracing steps. This includes numerically computing a location or position for a corresponding adjusted image surface or plane that corresponds to the required spherical dioptric power 1026 on which the image to be displayed will be mapped.

An exemplary ray-tracing methodology is described in steps 1104 to 1118 of FIG. 11 , at the end of which the output color of each pixel of pixel display 108 is known so as to virtually reproduce the light field emanating from an input image 1020 positioned at the adjusted image plane. In FIG. 11 , these steps are illustrated as a loop over each pixel in pixel display 108, so that each of steps 1104 to 1116 describes the computations done for each individual pixel. However, in some embodiments, these computations need not be executed sequentially, but rather, steps 1104 to 1116 may be executed in parallel for each pixel or a subset of pixels at the same time. Indeed, as will be discussed below, this exemplary method is well suited to vectorization and implementation on highly parallel processing architectures such as GPUs.

As illustrated in the schematic diagrams of FIG. 14A, in step 1104, for a given pixel 1402 in pixel display 108, a trial vector 1404 is first generated from the pixel's position to the pupil center position 1022 of the user's pupil 1405. As mentioned above, pupil center position 1022 may be acquired via eye/pupil tracker 110. In addition, a corresponding pupil plane 1406, which may be a flat 2D surface or plane in 3D space centered on pupil center position 1022, may be defined here as well. As shown in FIG. 14A, trial vector 1404, by construction, necessarily has to go through a corresponding optical unit of LFSL 106.

Once trial vector 1404 has been computed, in step 1106, a new Ray vector 1412 will be similarly generated from a center location/position 1410 of the corresponding optical unit comprising intersection point 1408 of LFSL 106 and pointing to pixel 1402. In this exemplary embodiment, step 1106 is detailed in the sub-steps 1206 to 1212 shown in FIG. 12 .

Thus, in sub-step 1206, the location of intersection point 1408 of vector 1404 with the LFSL 106 is calculated as illustrated in FIG. 14A. In sub-step 1208, the coordinates of the center location 1410 of the optical element or unit of LFSL 106 closest to intersection point 1408 are computed. One way to efficiently compute this location is provided in related U.S. Pat. No. 10,394, 322, the contents of which are incorporated herein by reference.

Once the position of the center 1410 of the optical element of LFSL 106 is known, in step 1210, as mentioned above, a normalized unit ray vector is generated from normalizing a ray vector 1412 originating from center position 1410 of LSFL 106 and extending to pixel 1402. This unit ray vector thus approximates the direction of the light field emanating from pixel 1402 through the center 1410 of this particular LFSL element, for instance, when considering a parallax barrier aperture or lenslet array (i.e. where the path of light travelling through the center of a given lenslet is not deviated by this lenslet). Further computation may be required when addressing more complex light shaping elements, as will be appreciated by the skilled artisan.

The orientation of ray vector 1412 will be used to find the portion of input image 1020 on the adjusted image plane, and thus the associated color, represented by pixel 1402. But first, in step 1212, ray vector 1412 is projected backwards (dotted line 1414 on FIG. 14B) to intersect with pupil plane 1406 at location 1416.

Going back to FIG. 11 , at step 1108, method 1100 verifies that intersection point 1416 with pupil plane 1406 of projected ray vector 1414 is still located within user pupil entrance 1405 (i.e. that the user can still “see” it). Thus, once intersection point 1416 shown in FIG. 14B of projected ray vector 1414 with the pupil plane 1406 is known, the distance between the pupil center 1022 and intersection point 1416 within pupil plane 1406 may be calculated to determine if the deviation is acceptable, for example by using pre-determined pupil size 1424 and verifying how far the projected ray vector intersection 1416 is from pupil center 1022 within pupil plane 1406.

If this deviation is deemed to be too large (i.e. light emanating from pixel 1402 channeled through optical element center 1410 is not perceived by pupil 1405), then in step 1110, method 1100 flags pixel 1402 as unnecessary and to simply be turned off or to render a black color.

In other embodiments, step 1108 may be modified so that instead of having a binary choice between the ray vector hitting the pupil or not, one or more smooth interpolation function (i.e. linear interpolation, Hermite interpolation or similar) may be used to quantify how far or how close intersection point 1416 is to pupil center 1022 within pupil plane 1406 and outputs a corresponding continuous value between 1 or 0. For example, the assigned value is equal to 1 substantially close to pupil center 1022 and gradually change to 0 as intersection point 1416 substantially approaches the pupil edges or beyond. In this case, the branch containing step 1110 may be ignored completely and step 1108 goes directly to step 1112. Then, at the end of step 1114, which will be discussed below, the pixel color value computed therein for pixel 1402 will be modified to be somewhere between the full color value identified therein or black, depending on the value of the interpolation function used at step 1108 (1 or 0).

In yet other embodiments, pixels found to illuminate a designated area around the pupil may still be rendered, for example, to produce a buffer zone to accommodate small movements in pupil location, for example, or again, to address potential inaccuracies, misalignments or to create a better user experience.

In the case where ray vector 1414 is within the pupil entrance (or if an interpolation function is used as discussed above), at step 1112, a corresponding image portion of input image 1020 located on the adjusted image plane and its corresponding color value are identified. As discussed above, two different but equivalent adjusted image planes may be used: a virtual image plane 1502 as shown schematically in FIG. 15 (i.e. positioned behind pixel display 108), or the retinal plane 1702 as shown in FIGS. 17A to 17D (i.e behind user pupil 1405). Correspondingly, each variation of steps 1002 and 1112 are listed in the flow diagrams of FIGS. 13 and 16 , respectively.

Thus, step 1102 is illustrated in FIG. 13 , in accordance with one embodiment. In the case where virtual image planes are selected for ray-tracing in sub-step 1301, then as mentioned above ray-tracing is done using a virtual image plane 1502 as illustrated in FIG. 15 . At sub-step 1302, the location or distance of virtual image plane 1502 from pupil plane 1406 may be computed as a function of spherical dioptric power 1026 (i.e. using the thin lens formula or other) and/or minimum reading distance 1028 (and/or related parameters). Then, at step 1304, input image 1020 is mapped onto virtual image plane 1502 so that its size of also scaled so to ensure that the perceived light field image correctly fills pixel display 108 when viewed by the distant user. An example is shown in FIG. 15 wherein virtual image plane 1502 is shown, as an example only, being located at a distance 1028 (in the z direction or depth) from pupil plane 1406, and the size of image 1020 is increased to avoid having the image as perceived by the user appear smaller than the display's size.

Continuing for the case of ray-tracing on virtual image plane 1502, the correct image portion is then identified in step 1112. So, as illustrated in FIG. 16 , in this case sub-step 1601 leads to sub-step 1602. As schematically illustrated in FIG. 15 , in sub-step 1602, ray vector 1412 is projected towards virtual image plane 1502 (shown as vector 1504 in FIG. 15 ) to find the position of the intersection point 1506. After this, the portion of image 1020 (and its associated colour channel) corresponding to intersection point 1506 on virtual image plane 1502 is identified.

As mentioned above, it may also be possible to use a retinal plane, herein defined as a 2D plane or surface substantially located at the location of the user's retina, as the adjusted image plane instead of virtual plane 1502.

Thus, in this case, illustrated schematically in FIGS. 17A to 17D, the adjusted image portion associated with a given pixel/subpixel is computed (mapped) on retinal plane 1702 instead of virtual image plane 1502 considered in the above example, again in order to provide the user with a designated image perception adjustment.

The skilled artisan will understand that a retinal plane may be defined in various ways. The exemplary embodiment described herein, the retinal plane 1702 is defined as a 2D plane located at a distance inside the eye equal to eye depth 1034 from the location of pupil center location 1022. It may also be taken to be parallel to pupil plane 1406, as illustrated in FIG. 17A to 17D, although it is not required to be so. It is meant to be an approximate location corresponding to the user's real retina.

Thus, in the case where ray-tracing is done on retina image plane 1702, in step 1102 as shown in FIG. 13 , sub-step 1301 leads to sub-step 1308, where a projected image center position on retinal plane 1702 is calculated. To do so, as illustrated in FIG. 17C, a vector 1704 is drawn originating from the center 1706 of pixel display 108 and passing through pupil center 1022. Vector 1704 is further projected beyond pupil plane 1405 onto retinal plane 1702, and the associated intersection point 1708 gives the location of the corresponding image center on retinal plane 1702. Once image center 1708 is known, in sub-step 1310 one can scale image 1020 to the x/y retina image size 1710, as illustrated schematically in FIG. 17D. In some embodiments, the required scaling may be computed by calculating the magnification of an individual pixel on retinal plane 1702, for example, which may be approximately equal to the x or y dimension of an individual pixel multiplied by the eye depth 1034 and divided by the absolute value of the distance to the eye (i.e. thus giving the magnification of the image portion created a pixel once focused by the eye lens on retinal plane 1702). An exemplary scaled inverted image 1712 on retinal plane 1702 is shown in FIG. 17D. Similarly, for comparison purposes, the input image 1020 may also normalized by the image x/y dimensions to produce a corresponding normalized input image 1714 having a width and height between −0.5 to 0.5 units, which may be compared to inverted scaled image 1712 which may also be similarly normalized.

In addition to sub-steps 1308 and 1310, sub-step 1312 is done independently to determine a location of a focal plane as produced by the user's eye for a given input value of spherical dioptric power 1026 (or minimum reading distance 1028). Thus, eye focal plane 1716 shown in FIGS. 17A and 17B is defined as the location where any light ray originating from optical unit center location 1410 would be focused by the user's eye. For a user with perfect vision, focal plane 1716 would be located at the same location as retinal plane 1702, but in the example shown in FIGS. 17A and 17B, as an example only, focal plane 1716 is located behind retinal plane 1702, which would be expected for a user with some form of farsightedness. The position of focal plane 1716 may be derived from the user's minimum reading distance 1028 or spherical dioptric power 1026, for example, by deriving therefrom the corresponding focal length of the user's eye. Other manually input or computationally or dynamically adjustable means may also or alternatively be considered to quantify this parameter.

Going back to FIG. 16 , step 1112 in the case of ray-tracing on retinal plane 1702 has sub-step 1601 leading to sub-step 1608, illustrated schematically in FIG. 17A, where a vector 1718 is drawn from optical unit center 1410 to pupil center 1022. Then, in sub-step 1610, vector 1718 is projected further behind pupil plane 1406 onto eye focal plane 1716 where intersection point 1720 is identified.

The skilled artisan will note that any light ray originating from optical unit center 1410, no matter its orientation, will also be focused onto intersection point 1720, to a first approximation. Therefore, in some embodiments, the location 1722 on retinal plane 1702 onto which light entering the pupil at intersection point 1416 will converge may be approximated, at sub-step 1612, by drawing a straight line between intersection point 1416 where projected ray vector 1414 hits pupil plane 1406 and focal point 1720 on focal plane 1716, as illustrated in FIG. 17B. The intersection of this line with retinal plane 1702 (retina image point 1722) is thus the location on the user's retina corresponding to the image portion that will be reproduced by corresponding pixel 1402 as perceived by the user. Therefore, at sub-step 1614, by comparing the relative position of retinal point 1722 with the overall position of the projected image on the retinal plane 1702, the relevant adjusted image portion associated with pixel 1402 may be computed.

For example, the image portion position 1724 relative to retina image center position 1726 in the scaled coordinates (scaled input image 1714) corresponds to the inverse (because the image on the retina is inverted) scaled coordinates of retina image point 1722 with respect to retina image center 1708, as shown in FIG. 17D. Thus, the associated color with image portion position 1724 may be therefrom extracted and associated with pixel 1402.

Once step 1112 is finished, in step 1114, pixel 1409 is flagged as having the color value associated with the portion of image corresponding to intersection point 1506 in the case of ray-tracing on virtual image plane 1502 (as shown in FIG. 15 ) or in the case of ray-tracing on retinal plane 1702, to the image portion corresponding to intersection 1725 as shown in FIG. 17D.

At step 1116, a check is made to see if every pixel in pixel display 108 has been ray-traced. If not then method 1100 chooses another pixel 1402 and goes back to step 1104; if so, then the output color of all pixels has been determined and these are finally rendered in step 1118 by pixel display 108 to be viewed by the user, therefore presenting a light field corrected image. In the case of a single static image, the method may stop here. In yet other embodiments, pixels found to illuminate a designated area around the pupil may still be rendered, for example, to produce a buffer zone to accommodate small movements in pupil location, for example, or again, to address potential inaccuracies or misalignments.

As will be appreciated by the skilled artisan, selection of the adjusted image plane onto which to map the input image in order to adjust a user perception of this input image allows for different ray tracing approaches to solving a similar challenge, that is of creating an adjusted image using the light field display that can provide an adjusted user perception, such as addressing a user's reduce visual acuity. While mapping the input image to a virtual image plane set at a designated minimum (or maximum) comfortable viewing distance can provide one solution, the alternate solution may allow accommodation of different or possibly more extreme visual aberrations. For example, where a virtual image is ideally pushed to infinity (or effectively so), computation of an infinite distance becomes problematic. However, by designating the adjusted image plane as the retinal plane, the illustrative process steps of FIGS. 16A and 16B can accommodate the formation of a virtual image effectively set at infinity without invoking such computational challenges. Likewise, while first order aberrations are illustratively described with reference to FIG. 11 , higher order or other optical anomalies may be considered within the present context, whereby a desired retinal image is mapped out and traced while accounting for the user's optical aberration(s) so to compute adjusted pixel data to be rendered in producing that image. These and other such considerations should be readily apparent to the skilled artisan.

While the computations involved in the above described ray-tracing algorithms (steps 1104 to 1116 of FIG. 11 ) may be done on general CPUs, it may be advantageous to use highly parallel programming schemes to speed up such computations. While in some embodiments, standard parallel programming libraries such as Message Passing Interface (MPI) or OPENMP may be used to accelerate the light field rendering via a general-purpose CPU, the light field computations described above are especially tailored to take advantage of graphical processing units (GPU), which are specifically tailored for massively parallel computations. Indeed, modern GPU chips are characterized by the very large number of processing cores, and an instruction set that is commonly optimized for graphics. In typical use, each core is dedicated to a small neighborhood of pixel values within an image, e.g., to perform processing that applies a visual effect, such as shading, fog, affine transformation, etc. GPUs are usually also optimized to accelerate exchange of image data between such processing cores and associated memory, such as RGB frame buffers. Furthermore, smartphones are increasingly being equipped with powerful GPUs to speed the rendering of complex screen displays, e.g., for gaming, video, and other image-intensive applications. Several programming frameworks and languages tailored for programming on GPUs include, but are not limited to, CUDA, OpenCL, OpenGL Shader Language (GLSL), High-Level Shader Language (HLSL) or similar. However, using GPUs efficiently may be challenging and thus require creative steps to leverage their capabilities, as will be discussed below.

In some embodiments, additional efficiencies may be leveraged on the GPU by storing the image data, for example image 1020, in the GPU's texture memory. Texture memory is cached on chip and in some situations is operable to provide higher effective bandwidth by reducing memory requests to off-chip DRAM. Specifically, texture caches are designed for graphics applications where memory access patterns exhibit a great deal of spatial locality, which is the case of the steps 1104 to 1116 of FIG. 11 . For example, in method 1100, image 1020 may be stored inside the texture memory of the GPU, which then greatly improves the retrieval speed during step 1112 (including either one of the ray-tracing variants presented in FIGS. 13 and 16 ) where the color channel associated with the portion of input image 1020 is determined.

Moreover, while method 1100 presented above (and its associated variations) was discussed and illustrated as having each plane (i.e. virtual image plane 1502, pixel display 108, LSFL 106, pupil plane 1406, retinal plane 1702 or eye lens focal plane 1716) as being parallel with each other, this was only done as an example for clarity and to better describe the methodology associated therewith. Indeed, method 1100 as discussed may equally be applied to account for changes in the relative orientation between any one of those planes.

For example, and with reference to FIGS. 18 and 19 , and in accordance with one exemplary embodiment, ray-tracing with non-parallel planes will now be discussed. In some embodiments, and as illustrated in FIG. 18 , cases may be considered wherein the user is viewing the light field display at an angle. In this specific example, the ray-tracing method can therefore account for a change in orientation of the pupil plane 1406 with respect to the pixel display 108 and optical layer 106. In this example, other planes such as virtual image plane 1502, and retinal plane 1702 and focal plane 1716 may be taken to be parallel to pupil plane 1406. The relative difference in orientation between the two sets of planes is illustrated by using vector 1802 which is the normal vector to the plane of corresponding optical layer 106, and vector 1804 which is the normal vector to pupil plane 1406. The relative orientation between the two normal vectors is illustrated in FIG. 19 , using polar and azimuthal angles.

The general orientation of pupil plane 1406 may be parametrized, for example, by using the 3D location of pupil center 1022 and a corresponding normal vector 1804. Normal vector 1804 may be taken to be, in some embodiments, equal to the gaze direction as measured by a gaze tracking system or similar, as will be discussed below.

Once the relative position and orientation of pupil plane 1406 is determined, the relative position/orientation of all remaining planes (parallel or non-parallel) may be determined and parametrized accordingly. Planes that are parallel share the same normal vector. From there, the method 1100 and its variants described above may be applied by finding the intersection point between an arbitrary vector and an arbitrarily oriented plane, as is done for example at steps 1206, 1212, 1602, 1610, 1612 for example.

In the illustrated example of FIG. 18 , the position of virtual image plane 1502 may be computed using spherical dioptric power 1026 (and/or minimum reading distance 1028 and/or related parameters) but from the position of pupil plane 1406 and along the direction vector 1804.

To extract normal vector 1804 of pupil plane 1406, the eye tracking methods and systems described above may be used or modified to further provide a measure of the eye's gaze direction (e.g. gaze tracking). As discussed above, there are many known eye tracking methods in the art, some of which may also be used for gaze-tracking. For example, this includes Near-IR glint reflection methods and systems or methods purely based on machine vision methods. Hence, in some embodiments, pupil plane 1406 may be re-parametrized using an updated 3D location of pupil center 1022 and an updated normal vector 1804 at each eye tracking cycle. In other embodiments, a hybrid gaze tracking/pupil tracking system or method may be used wherein gaze direction (e.g. normal vector 1804) is provided at a different interval than pupil center location 1022. For example, in some embodiments, for one or more cycles, only the 3D pupil center location 1022 may be measured and an old gaze direction vector may be re-used or manually updated. In some embodiments, an eye model or similar may be constructed to map a change in measured pupil center location 1022 to a change in the gaze direction vector without relying on the full capabilities of the gaze tracking system or method. Such a map may be based on one or more previous gaze tracking measurements. In any case, by measuring/determining the 3D pupil center location 1022 and normal vector 1804, the pupil plane may be parametrized accordingly.

Note that in FIG. 18 , pixel display 108 and optical layer 106 are shown as being parallel for simplicity, but other embodiments may envision optical layer 106 to be non-parallel to display 108 as well. This doesn't change the general scope of the present discussion, as long as the relative angle between them is known. For example, such an angle may be pre-determined during manufacturing or measured in real-time using one or more sensors (for example in the case where optical layer 106 may be mobile). Similarly, other planes like for example retinal plane 1702 may also be made to be non-parallel to the pupil plane, depending on the user's eye geometry.

With reference to FIGS. 20A to 22D and in accordance with one embodiment, a modified embodiment of method 1100 operable to render multiple images or image portions on multiple adjusted distinct image planes simultaneously via an array of light field shaping elements, or light field shaping layer (LFSL) thereof, will now be discussed. Thus, step 2002 shown in FIG. 20A is meant to replace step 1102 in method 1100, while step 2012 of FIG. 20B replaces step 1112. This is because the previously above-described steps 1102 and 1112 were directed to correcting a single image by directly or indirectly modifying the location of the virtual image plane and/or eye focal plane. In contrast, the below-described embodiments of steps 2002 and 2012 are directed to a light field display which is generally operable to display multiple image planes at different locations/depths/aberrations simultaneously. Thus, in step 2002, at sub-step 2004, the method continues to sub-step 2102 (illustrated in FIG. 21A) if ray-tracing is done using the virtual image plane or if ray-tracing is done using the retinal plane then sub-step 2300 is used (illustrated in FIG. 23A). Similarly, as shown in FIG. 20B, from sub-step 2014, the method proceeds to sub-step 2104 in the case of ray-tracing to virtual image plane (illustrated in FIG. 21B) or to sub-step 2301 (illustrated in FIG. 23B if ray-tracing is done using the retinal plane).

Unlike known stereoscopic effects, the methods as herein described may be implemented to generate varying depth perceptions within a same eye, that is, allowing for the monoscopic viewing of an input to exhibit multiple distinct image perception adjustments (i.e. multiple juxtaposed and/or overlapping depths, enhancements or like optical adjustments, compensations, etc.). For example, in some embodiments, distinct image planes may be juxtaposed such that different sides or quadrants of an image, for example, may be perceived at different depths. In such embodiments, a different effective vision correction parameter (e.g. diopter), or depth, may be applied, to each portion or quadrant. While this approach may result in some distortions or artefacts at the edges of the areas or quadrants, depending on the image data to be rendered along these edges, such artefacts may be negligible if at all perceivable. In other embodiments, however, different image portions may be at least partially superimposed such that portions at different depths, when viewed from particular perspectives, may indeed appear to overlap. This enables a user to focus on each plane individually, thus creating a 2.5D effect. Thus, a portion of an image may mask or obscure a portion of another image located behind it depending on the location of the user's pupil (e.g. on an image plane perceived to be located at an increased distance from the display than the one of the first image portion). Other effects may include parallax motion between each image plane when the user moves.

As mentioned above, steps 2102 and 2104 of FIGS. 21A and 21B are directed towards ray-tracing on one or more virtual image plane only, while steps 2300 and 2301 of FIGS. 23A and 23B are directed towards ray-tracing on the retinal plane.

For example, to account for multiple distinct image planes, input image 1020 of input variables 1004 may also include, in addition to pixel data, variable dioptric powers or perceptive “depth” information or parameters. Thus, any image or image portion may have a respective depth indicator. Thus, at sub-step 2106, a set of multiple virtual image planes may be defined, at sub-step 2108, which includes deriving their respective (virtual) location, similarly to sub-step 1302. On these planes, images or image portions may be present. Areas around these images may be defined as transparent or see-through, meaning that a user would be able to view through that virtual image plane and see, for example, images or image portions located behind it. At sub-step 2108, any image or image portion on each of these virtual image planes may be optionally scaled to fit the display, similarly as described for sub-step 1304 for a single image plane.

In the previous example shown in FIG. 15 , a single virtual image plane 1502, showing an exemplary input image 1020 comprising two circles, was used. In contrast, FIGS. 22A to 22D show an example wherein each circle is located on its own virtual image plane (e.g. original virtual plane 1502 with new virtual image plane 2202). The skilled technician will understand that two planes are shown here only as an example and that the method steps described herein apply equally well to any number of virtual planes. The only effect of having more planes is a larger computational load.

Going back to 21B, in step 2104, an iteration is done over the set of virtual image planes to compute which image portion from which virtual image plane is seen by the user. Thus, at sub-step 2110 a virtual image plane is selected, starting from the plane located closest to the user. Then step 1602 proceeds as described previously for that selected virtual plane. At sub-step 2112 the corresponding color channel of the intersection point identified at step 1602 is sampled. Then at sub-step 2114, a check is made to see if the color channel is transparent. If this is not the case, then the sampled color channel is sent to step 1114 of FIG. 11 , which was already described and where the color channel is rendered by the pixel/subpixel. An example of this is illustrated in FIGS. 22A and 22B, wherein a user is located so that a ray vector 2204 computed passing through optical element center location 1410 and pixel/subpixel 1402 intersects virtual image plane 1502 at location 1506. Since in this example this location is non-transparent, this is the color channel that will be assigned to the pixel/subpixel. However, as this example shows, this masks or hides parts of the image located on virtual image plane 2202. Thus, an example of the image perceived by the user is shown in FIG. 22B.

Going back to FIG. 21B, at sub-step 2114 if the color channel is transparent, then another check is made at sub-step 2116 to see if all virtual image planes have been iterated upon. If this is the case, then that means that no image or image portion is seen by the user and at sub-step 2118, for example, the color channel is set to black (or any other background colour), before proceeding to step 1114 of FIG. 11 . If however at least one more virtual image plane is present, then the method goes back to step 2110 and selects that next virtual image plane and repeats sub-steps 1602, 2112 and 2114. An example of this is illustrated in FIG. 22C, wherein a user is located so that a distinct ray vector 2204 computed passing through optical element center 2206 of LFSL 106 and pixel/subpixel 2208 of pixel display 108 first intersects at location 2210 of virtual image plane 1502. Since, in this example, this location is defined to be transparent (i.e. not on the circle), the method checks for additional virtual image planes (here plane 2202) and extends vector 2204 so as to compute intersection point 2212, which is now non-transparent (i.e. on the circle), and thus the corresponding color channel is selected. An example of the image perceived by the user is shown in FIG. 22D.

Similarly, steps 2300 and 2301 of FIGS. 23A and 23B substantially mirrors steps 2102 and 2104, respectively, described in FIGS. 21A and 21B, but are herein applied to be used with two or more eye focal planes (e.g. for ray-tracing the image on retinal image plane 1702). Thus, we see that the method iterates, after sub-step 1608, over all designated image planes, each corresponding to a different eye focal plane, starting from the plane corresponding to an image located closest to the user. Thus, a new eye focal plane is selected at sub-step 2306, which is used for sub-steps 1610 to 1614 already described above. Once the corresponding image portion is determined at sub-step 1614, at sub-step 2308, the corresponding pixel/subpixel color channel is sampled. Then at sub-step 2310, if the color channel is non-transparent, then the method continues to step 1114 of FIG. 11 , wherein the pixel/subpixel is assigned that color channel. However, if the image portion is transparent, then the method iterates to the eye focal plane corresponding to the next designated image plane. Before this is done, the method checks at sub-step 2312 if all the eye focal planes have been iterated upon. If this is the case, then no image portion will be selected and at sub-step 2314 the color channel is set to black, for example, before exiting to step 1114. If other eye focal planes are still available, then the method goes back to sub-step 2306 to select the next eye focal plane and the method iterates once more.

In some embodiments, as mentioned above, steps 2102, 2104 or 2300 and 2301 for multiple designated image planes of FIGS. 21A-B (on virtual planes) or FIGS. 23A-B (retinal plane) may be used to implement a phoropter/refractor device to do subjective visual acuity evaluations. For example, as illustrated in FIGS. 24A and 24B, and similarly to FIG. 9 , different optotypes (e.g. letters, symbols, etc.) may be displayed simultaneously but at different perceived depths, to simulate the effect of adding a refractive optical component (e.g. change in focus/optical power). In FIG. 24A, two images of the same optotype (e.g. letter E) are displayed, each on their own designated image plane (e.g. here illustrated as virtual image planes as an example only). In this example, image 2402 is located on designated image plane 2404 while image 2406 is located on designated image plane 2408, which is located further away. Optionally, as illustrated herein, the size of the image may be increased with increased depth so that all images displayed are perceived to be of a similar relative size by the user. In FIG. 24B, we see an example of the perception of both images as perceived by a user with reduced visual acuity (e.g. myopia), for example, wherein the image closest to the user is seen to be clearer. Thus, a user could be presented with multiple images (e.g. 2 side-by-side, 4, 6 or 9 in a square array, etc.) and indicate which image is clearer and/or and most comfortable to view. An eye prescription may then be derived from this information.

As mentioned above, the systems and methods discussed herein may also be adapted to correct or compensate for the presence of HOAs. This is generally done by digitally simulating the presence of a “virtual” optical element specifically tailored or constructed so as to reproduce the effect of said HOAs. Such a “virtual” optical element will be referred to herein as a phase element, since its effect on an incoming wave front is to introduce a phase shift or phase profile corresponding to one or more known HOAs. Thus, by considering the effects of such a phase element on each incoming ray generated by a light field display, it is possible to generate an image that compensates for the user's HOAs and thus allows the user to clearly perceive the image. Notably, although the focus of the methods described below is one treating for HOAs, the methods using the phase element may generally be applied, without limitation, to simulate other kinds of optical aberrations, such as astigmatism, as will be made clear below.

The use of a virtual phase element is schematically illustrated in FIG. 25 , wherein an incoming source of light herein pictured as a plane wave front 2502 is refracted by a phase element 2500, which is tailored or constructed so as to generate an aberrated output wave front 2504. Since a wave front is defined as asurface of common phase, a corresponding set of light rays may be defined so as to be perpendicular or normal to the wave front surface (incoming rays 2506 and refracted rays 2508). The difference between aberrated output wave front 2504 and input wave front 2502 is in the phase shift caused by one or more HOAs (or aberrations in general) of phase element 2500. Thus, by constructing or defining phase element 2500 so as to reproduce the required phase shift, it is possible to compute therefrom the change in orientation of input rays 2506 to output rays 2508 and construct a light field image using a light field display that compensates for these HOAs (or aberrations in general).

Notably, phase element 2500 is herein illustrated as being parallel to the (x, y) axis as an example only. More general relative orientations between phase element 2500 and for example light field display 104 or pupil plane 1406 (not shown here) may be assumed, without limitation, as will be noted further below.

As schematically illustrated in FIGS. 26A and 26B, in some embodiments, phase element 2500 may be used to model or digitally replicate HOAs (or aberrations in general) in the user's eye, in which case phase element 2500 may be virtually positioned at the user's pupil, as illustrated in FIG. 26A. In other embodiments, it may be desirable to model HOAs (or wavefront modifications) of a real physical optical element located in front of the eye, and thus phase element 2500 may be positioned somewhere in-between the eye and the light field display, as shown in FIG. 26B. In some embodiments, multiple phase elements 2500 may be deployed simultaneously.

Phase element 2500 may be defined or constructed in different ways and may have different geometries or optical characteristics, so long as it correctly models or reproduces the corresponding aberrations of interest. Some non-limiting examples are given below.

For example, as seen in FIG. 27A, in some embodiments, phase element 2500 may comprise a first incident surface 2700 (facing the light field display), a second surface 2700 (facing the user). However, different combinations of surface profiles may be considered, without limitation. For example, in some embodiments, as in the illustrated exemplary embodiment of FIG. 27B, first incident surface 2700 may be defined to be flat with a normal vector orientated in the z direction, and only the second surface 2702 having a variable shape profile. Other embodiments may include both surfaces having a variable shape or profile, as illustrated in FIG. 27C.

In addition, phase element 2500 is not limited to having a rectangular shape, as was illustrated in FIGS. 26A-C and 27A only as an example. Different considerations may be taken into account to determine the shape, size or extension of phase element 2500 in the (x, y) and z dimensions. In some embodiments, for example when Zernike polynomials are used to generate the shape of phase element 2500, a disk shape, as in FIG. 27D, may be used, so as to better align with the pupil shape and the azimuthal and radial dependent forms of the Zernike polynomials. Generally, the size in the (x, y) plane will be chosen so that all light rays projected from the light field display may be properly refracted. For example, in the case where phase element 2500 is placed at the pupil plane 1406, (discussed further below), the x, y coordinates may be chosen as being equal to a factor of the pupil radius, so as to ensure that all rays encountering the pupil are considered. For example, as illustrated in FIG. 27D, different in-plane shapes (herein shown as a disk-shape or circular profile) may be considered, without limitation, as required.

In some embodiments, phase element 2500 may comprise multiple layers (not shown), each layer having its own refractive index, thus comprising a multiplicity of internal surfaces in addition to the external surfaces 2700 and 2702 (not shown).

In the methods considered below, as an example only, the embodiment of phase element 2500 considered, will consist of a phase element having a circular profile (as in the diagram of FIG. 27D), but having a flat or plane-like first surface 2700 and a second surface 2702 having a variable shape, depth, height or profile defined via a height or thickness function h(x, y), as illustrated schematically in FIG. 28 . In this example, function h(x, y) may thus give the height or depth z-position of second surface 2702 as measured perpendicularly from the z-position of first surface 2700 for all (x, y) values. In this example, phase element 2500 is thus constructed by computing the height or thickness function h(x, y) (in conjunction with a refractive index n) that reproduces the required HOAs (or aberrations in general).

Different methods may be used to derive the correct profile of phase element 2500. In some embodiments, the required aberrational phase shift may be defined using Zernike Polynomials, for example. In other embodiments, the required phase shift may be derived from measurements or the like, for example a lens profile measurement or wave front measurement. The embodiments discussed below will use Zernike polynomials, as an example only, since they are commonly used as a complete set of functions to represent the effect of any optical element or aberration on an incoming wave front. However, other decomposition functions or polynomials may be used, for example and without limitation, Seidel series, Fourier series, Taylor polynomials, elliptical polynomials, etc.

In addition, in some embodiments, small beams may be assumed when sampling the phase element 2500, so that the divergence or convergence due to the presence of the phase element may be minimal. Again, the various embodiments described below will use such an approximation, also as an example only.

Now, with reference to FIGS. 29 to 35 , and in accordance with different exemplary embodiments, a method for correcting or compensating for HOAs (or aberrations in general) via a light field display, generally referred to using the numeral 2900, will now be discussed.

In some embodiments, method 2900 as illustrated in FIG. 29 is an extension of method 1100 discussed above, but now includes the effects of incorporating one or more phase elements 2500 to treat HOAs (or aberrations in general). While some steps are similar as those described above, other steps have been modified or added as required.

At step 2902, the image or optotype to be displayed is prepared. This is illustrated in FIG. 30 , where, at sub-step 3002, one or more virtual phase elements 2500 are constructed or defined. This includes, in part, defining the required geometry and optical properties of each phase element 2500 so as to reproduce the required aberrational phase shift, as was explained above. Notably, in some embodiments, multiple phase elements 2500 may be used simultaneously, as will be discussed below.

An example of sub-step 3002 is shown in FIG. 31A. There, an iteration over all the phase element(s) 2500 to be constructed or defined is made starting at sub-step 3101. For each phase element 2500, its geometry and refraction index are computed at sub-step 3102 and its orientation and location at sub-step 3104. At sub-step 3105 the method proceeds back to sub-step 3101 if more phase elements are to be used.

In one exemplary embodiment, and as illustrated in FIG. 31B, a given phase element 2500 is defined or constructed in sub-step 3102 to reproduce a phase shift as defined by a given Zernike polynomial which may be characterized by a corresponding set of Zernike coefficients 1036. In addition, in this embodiment, it is designed so as to comprise a first input surface 2700 (facing the light field display) that is flat and parallel to the (x, y) plane (as an example only) and second output surface 2702 (facing the user) which has a profile or shape that varies as a function of planar position (x, y). In addition, a refractive index 1038 (n) of phase element 2500 is also chosen accordingly. The exact value of refractive index 1038 is not meant to reproduce exactly the refractive index of a real optical element (external lens or user's eye), but it is chosen so that, in combination with the height function h(x, y), the phase element 2500 properly or correctly reproduces the chosen HOAs (or aberration in general) evaluated by comparison to diffraction patterns of wavefronts (for example such as those illustrated in FIGS. 32A and 32B).

Notably, if the height or thickness function h(x, y) is known analytically, corresponding analytical expressions for the normal vector functions of the second surface 2702 may also be used. In other embodiments, height function h(x, y) and refractive index/indices may be derived from experimental data or similar, and thus they can be used directly, and it may be required to use numerical derivatives or similar methods to compute the normal vectors. In addition, the data points can be fitted to smooth analytical functions to reduce the numerical errors in calculating the derivatives.

In some embodiments, the pre-computed or measured height or thickness function h(x, y) of phase element 2500 and the associated pre-computed normal vector coordinates may be stored digitally in an array structure or similar. For example, these values may be stored via a texture object or similar, for ease and speed of access. Thus, in some embodiments, different “resolutions” of the texture object may be used, so as to “coarse-grain” the geometry as required by considerations of speed of computation, as will be further discussed below.

In one exemplary embodiment of step 3102, illustrated in FIG. 31B, the HOAs are decomposed or described mathematically using Zernike polynomials (Z_(i) or Z_(k) ^(j)). These polynomials are used to represent or define deviations or shifts from a desired wave-front. Since a wave-front is a plane of uniform phase, the Zernike phase shift describing the effect of one or more optical aberrations may be written as, for example:

${\varphi_{Z}\left( {x,y} \right)} = {{{\mathbb{Z}}\left( {x,y} \right)} = {{\frac{2\pi}{\lambda_{0}}{\sum{C_{i}{Z_{i}\left( {x,y} \right)}}}} = {\frac{2\pi}{\lambda_{0}}{\sum{C_{kj}{{Z_{k}^{j}\left( {r,\theta} \right)}.}}}}}}$

In some embodiments, as was illustrated in FIGS. 28A and 28B, the spatial dependent height or thickness h(x, y) of phase element 2500, as defined from the location or position of the first interface or surface 2700 with respect to the second surface 2702, may be computed using a corresponding phase shift function, φ(x, y), which represents the effect of HOAs (or aberrations in general) on an incoming wave front, and at least one corresponding refractive index 1038 (n) of phase element 2500. In some embodiments, the phase shift may be constructed on a plane parallel to the lens principal plane using a multiplicity of parallel rays passing through the interior of phase element 2500, again as was illustrated in FIG. 25 . To achieve that, both phase element 2500 and the remnant distance to reach the phase-construction plane may be considered.

Thus, as mentioned above, at sub-step 3102 of FIG. 31A, the geometry and refraction index of one phase element 2500 is generated from the input set of Zernike coefficients 1036. For example, in the exemplary embodiment of FIGS. 28A and 28B, a thin element assumption or approximation may be used, so that the phase-change or phase-shift as a function of (x, y) coordinates generated by phase element 2500 may be written as:

${\varphi\left( {x,y} \right)} = {\frac{2\pi}{\lambda_{0}}\left\lbrack {{{nh}\left( {x,y} \right)} + {n_{0}\left( {h_{\max} - {h\left( {x,y} \right)}} \right)}} \right\rbrack}$

which, excluding the constant phase term, gives a height or thickness as a function of the (x, y) coordinates of:

${h\left( {x,y} \right)} = {\frac{\lambda_{0}}{2\pi}\frac{\varphi\left( {x,y} \right)}{n - n_{0}}}$

where n₀ is the medium refractive index and n is the refractive index 1038. Assuming phase element 2500 is immersed in air (i.e. n₀=1), then the (height or thickness) becomes:

${h\left( {x,y} \right)} = {\frac{\lambda_{0}}{2\pi}{\frac{\varphi\left( {x,y} \right)}{n - 1}.}}$

Thus, for a given phase-shift function φ(x, y), the equation above may be used to derive the corresponding height or thickness function h(x, y) at any point (x, y) of phase element 2500.

In addition, the refraction properties of phase element 2500 may be derived in different ways, as required. For example, and in accordance with one embodiment, the continuity of the tangential wave-vector component may be used to derive ray refraction equations. The position (x_(T), y_(T)) of the transmitted ray through a surface is the same as the incident ray (x_(I), y_(I)), thus: (x _(T) ,y _(T))=(x _(I) ,y _(I))

A normal of the surface vector (N _(S)) and a corresponding incident ray wave-vector (K _(I)) define the plane of incidence. This plane can be defined by the normal that is calculated by the cross product:

${\hat{N}}_{I} = {\frac{{\overset{\_}{K}}_{I} \times {\overset{\_}{N}}_{S}}{{❘{\overset{\_}{K}}_{I}❘}{❘{\overset{\_}{N}}_{S}❘}}.}$

The angle of incidence (θ_(I)) between the surface normal and the ray wave-vector, is given by:

${{\cos\left( \theta_{I} \right)} = \frac{{\overset{\_}{K}}_{I} \cdot {\overset{\_}{N}}_{S}}{{❘{\overset{\_}{K}}_{I}❘}{❘{\overset{\_}{N}}_{S}❘}}},$

So, in the plane of incidence, Snell's law can be used to calculate the transmitted ray angle (θ_(T)): n _(I) sin(θ_(I))=n _(T) sin(θ_(T))

where n_(I), n_(T), are the refractive indices in the incidence and transmission media, respectively.

In addition, the incident wave vector can be decomposed in terms of the surface normal vector and a surface tangential vector (N _(t)):

${a_{tI}{\overset{\_}{N}}_{t}} = {{\overset{\_}{K}}_{I} - {{\cos\left( \theta_{I} \right)}{❘{\overset{\_}{K}}_{I}❘}{\frac{{\overset{\_}{N}}_{S}}{❘{\overset{\_}{N}}_{S}❘}.}}}$

The transmitted ray direction-vector (K _(T)) may be calculated by considering that the continuity of the tangential component of the wave vector and the change in its magnitude. Again, representing the transmitted component as normal and tangential components, in the plane of incidence:

${\overset{\_}{K}}_{T} = {{\frac{{\overset{\_}{N}}_{S}}{❘{\overset{\_}{N}}_{S}❘}{\cos\left( \theta_{T} \right)}{❘{\overset{\_}{K}}_{I}❘}\frac{n_{T}}{n_{I}}} + {\overset{\_}{K}}_{I} - {{\cos\left( \theta_{I} \right)}{❘{\overset{\_}{K}}_{I}❘}{\frac{{\overset{\_}{N}}_{S}}{❘{\overset{\_}{N}}_{S}❘}.}}}$

In some embodiments, as in this example wherein the function h(x, y) is available analytically, corresponding analytical expressions for the normal unit vectors on the surface at every location h(x, y) may be derived using the gradient of the surface. Assuming the height or thickness is in the z direction only, for a function of the form: S(x,y,z)=0, the normal unit vectors to the surface are given by:

${\hat{n} = {{{\frac{\left\langle {S_{x},S_{y},S_{z}} \right\rangle}{\left. {❘\left\langle {S_{x},S_{y},S_{z}} \right.❘} \right\rangle}{for}\left( {n - n_{0}} \right)z} - {\frac{\lambda_{0}}{2\pi}{\varphi\left( {x,y} \right)}}} = 0}},{{{then}:\hat{n}} = {\frac{{\overset{\_}{N}}_{S}}{❘{\overset{\_}{N}}_{S}❘} = \frac{\left\lbrack {{- \frac{\lambda_{0}{\partial\varphi}}{2\pi{\partial x^{\prime}}}} - {\frac{\lambda_{0}{\partial\varphi}}{2\pi{\partial y^{\prime}}}\left( {n - n_{0}} \right)}} \right\rbrack}{\sqrt{\left( \frac{\lambda_{0}{\partial\varphi}}{2\pi{\partial x}} \right)^{2} + \left( \frac{\lambda_{0}{\partial\varphi}}{2\pi{\partial y}} \right)^{2} + \left( {n - n_{0}} \right)^{2}}}}},$

where n₀ is the medium refractive index and n is the phase element refractive index 1038. In air, this gives:

${\hat{n}\left( {x,y} \right)} = {\frac{{\overset{\_}{N}}_{S}}{❘{\overset{\_}{N}}_{S}❘} = {\frac{\left\lbrack {{- \frac{\lambda_{0}{\partial\varphi}}{2\pi{\partial x^{\prime}}}} - {\frac{\lambda_{0}{\partial\varphi}}{2\pi{\partial y^{\prime}}}\left( {n - 1} \right)}} \right\rbrack}{\sqrt{\left( \frac{\lambda_{0}{\partial\varphi}}{2\pi{\partial x}} \right)^{2} + \left( \frac{\lambda_{0}{\partial\varphi}}{2\pi{\partial y}} \right)^{2} + \left( {n - 1} \right)^{2}}}.}}$

In some embodiments, a variety of optimizations may be employed. As mentioned above, analytical forms of the derivatives may be used. In addition, a high value of refractive index 1038 may also be used to minimize the thickness of phase element 2500 (although precautions must be taken so as to avoid the total internal reflection regime).

In some embodiments, if the computed height or thickness function h(x, y) results in negative values and a non-thin phase element is considered, a constant height or thickness value may be subtracted before computing ray propagation inside the element. For example: h(x,y)=h(x,y)−min(h(x,y)).

where min(h(x, y)) is the minimum value of function h(x, y).

We present herein an example of how Zernike polynomials and their corresponding Zernike coefficients 1036 may be used to model HOAs.

In some embodiments, a thin lens assumption may be used for phase element 2500. For example, considering a thin element with a first flat incident surface 2700 is perpendicular to z-axis, then: S(x,y,z)=z−z ₀₌₀ {circumflex over (n)}=

0,0,1

and the second surface 2702, herein modelled as an ideal thin lensing surface may (with no approximation) may be given by:

${\phi\left( {x,y} \right)} = {{\frac{2\pi}{\lambda_{0}}\left\lbrack {f - \sqrt{x^{2} + y^{2} + f^{2}}} \right\rbrack} = {{\frac{2\pi}{\lambda_{0}}{f\left( {1 - \sqrt{1 + \frac{x^{2} + y^{2}}{f^{2}}}} \right)}} = {\frac{2\pi}{\lambda_{0}}\frac{1 - \sqrt{1 + {A^{2}\left( {x^{2} + y^{2}} \right)}}}{A}}}}$

where f is the focal length of phase element 2500;

which, for the paraxial limit gives:

$\begin{matrix} {{{\phi\left( {x,y} \right)} \approx {- \frac{2\pi}{\lambda_{0}}\frac{x^{2} + y^{2}}{2f}}} = {- \frac{2\pi}{\lambda_{0}}A\frac{x^{2} + y^{2}}{2}}} \\ {{S\left( {x,y,z} \right)} = {{{\left( {n - 1} \right)z} - f + \sqrt{x^{2} + y^{2} + f^{2}}} = 0}} \\ {\hat{n} = {\frac{\left\langle {\frac{x}{\sqrt{x^{2} + y^{2} + f^{2}}},\frac{y}{\sqrt{x^{2} + y^{2} + f^{2}}},\left( {n - 1} \right)} \right\rangle}{\sqrt{\frac{x^{2} + y^{2}}{x^{2} + y^{2} + f^{2}} + \left( {n - 1} \right)^{2}}}.}} \end{matrix}$

In the case of a thin lens in the paraxial approximation, transfer matrix model may be used. In this example, the spatial coordinates do not change, on the other hand the wave-vector angle is given by:

${{\overset{\_}{\theta}}_{T} = {\left( {\frac{- x_{I}}{f} + \theta_{xI}} \right) + \left( {\frac{- y_{I}}{f} + \theta_{yI}} \right)}}{{{For}:{\overset{\_}{K}}_{I}} = {{{a_{ix}\hat{x}} + {a_{iy}\hat{y}} + {a_{iz}\hat{z}\theta_{Ix}}} = {{{asin}\left( \frac{a_{ix}}{\sqrt{a_{ix}^{2} + a_{iz}^{2}}} \right)} \approx \frac{a_{ix}}{a_{iz}}}}}$ The transmitted directional vector

$\begin{matrix} {{\overset{\_}{K}}_{T} = {{a_{Tx}\hat{x}} + {a_{Ty}\hat{y}} + {a_{Tz}\hat{z}}}} \\ {\theta_{Tx} = {{{{asin}\left( \frac{a_{Tx}}{\sqrt{a_{Tx}^{2} + a_{Tz}^{2}}} \right)} \approx \frac{a_{Tx}}{a_{Tz}}} = {{\theta_{xI} - \frac{x_{I}}{f}} = {\frac{a_{ix}}{a_{iz}} - \frac{x_{I}}{f}}}}} \\ {a_{Tx} = {\left( {\frac{a_{ix}}{a_{iz}} - \frac{x_{I}}{f}} \right)a_{Tz}}} \\ {a_{Ty} = {\left( {\frac{a_{iy}}{a_{iz}} - \frac{y_{I}}{f}} \right)a_{Tz}}} \\ {\sqrt{a_{Tx}^{2} + a_{Ty}^{2} + a_{Tz}^{2}} = {❘{\overset{\_}{K}}_{T}❘}} \\ {{{a_{Tz}\sqrt{\left( {\frac{a_{ix}}{a_{iz}} - \frac{x_{I}}{f}} \right)^{2} + \left( {\frac{a_{iy}}{a_{iz}} - \frac{y_{I}}{f}} \right)^{2} + 1}} = {❘{\overset{\_}{K}}_{T}❘}},{{so}{that}}} \\ {{a_{Tz} = \frac{❘{\overset{\_}{K}}_{T}❘}{\sqrt{\left( {\frac{a_{ix}}{a_{iz}} - \frac{x_{I}}{f}} \right)^{2} + \left( {\frac{a_{iy}}{a_{iz}} - \frac{y_{I}}{f}} \right)^{2} + 1}}};} \\ {{a_{Ty} = \frac{{❘{\overset{\_}{K}}_{T}❘}\left( {\frac{a_{iy}}{a_{iz}} - \frac{y_{I}}{f}} \right)}{\sqrt{\left( {\frac{a_{ix}}{a_{iz}} - \frac{x_{I}}{f}} \right)^{2} + \left( {\frac{a_{iy}}{a_{iz}} - \frac{y_{I}}{f}} \right)^{2} + 1}}};} \\ {a_{Tx} = {\frac{{❘{\overset{\_}{K}}_{T}❘}\left( {\frac{a_{ix}}{a_{iz}} - \frac{x_{I}}{f}} \right)}{\sqrt{\left( {\frac{a_{ix}}{a_{iz}} - \frac{x_{I}}{f}} \right)^{2} + \left( {\frac{a_{iy}}{a_{iz}} - \frac{y_{I}}{f}} \right)^{2} + 1}}.}} \end{matrix}$

More generally, as discussed above, any phase shift φ(x, y) may be represented via a Zernike polynomial in the form of:

${\varphi\left( {x,y} \right)} = {{{\mathbb{Z}}\left( {x,y} \right)} = {\frac{2\pi}{\lambda_{0}}{\sum{C_{i}{Z_{i}\left( {x,y} \right)}}}}}$

and where the normal vectors to the non-flat surface of the corresponding phase element, in air, are given by:

$\hat{n} = {\frac{\left\langle {{- {\sum{C_{i}\frac{\partial{Z_{i}\left( {x,y} \right)}}{\partial x}}}},{- {\sum{C_{i}\frac{\partial{Z_{i}\left( {x,y} \right)}}{\partial y}}}},\left( {n - 1} \right)} \right\rangle}{\sqrt{\left( {\sum{C_{i}\frac{\partial{Z_{i}\left( {x,y} \right)}}{\partial x}}} \right)^{2} + \left( {\sum{C_{i}\frac{\partial{Z_{i}\left( {x,y} \right)}}{\partial y}}} \right)^{2} + \left( {n - 1} \right)^{2}}}.}$

Thus, for a given set of Zernike coefficients and their corresponding Zernike function, the corresponding height or thickness function h(x, y) and a corresponding set of normal vector functions {circumflex over (n)}(x, y) for the second surface may be computed.

In some embodiments, the refractive index and shape of phase element 2500 may be obtained by trial and error to match, for example, a refracted beam produced by phase element 2500 with a diffraction profile or by profiling a real optical element. For example, FIGS. 32A and 32B shows a series of plots illustrating the different phase element intensity profiles for refracted rays through a phase element and corresponding profiles for diffracted phase fronts for different Zernike coefficients (numbered via the OSA/ANSI standard indices). These correspond to a phase element in air (n_(o)=1) having a refractive index n=100, further assuming that the phase front focal length of 100 mm and a distance between the phase element and the chosen image plane of 1 m.

More generally, in the table below is an exemplary list of Zernike functions and their partial derivatives for different HOAs. Therein the variable r is the polar radial coordinates normalized to the pupil radius (R_(P)) and θ is the polar angular coordinate and z is the corresponding contribution to the height or thickness h(x, y).

Aberration OSAANSI (i) Noll Index Radial Index (k) Azimuthal Index (j) Z $\frac{\partial Z}{\partial x}$ $\frac{\partial Z}{\partial y}$ Piston 0 1 0   0 1 0 0 Y-Tilt 1 3 1 −1 2r sin(θ) = 0 2/R_(p) 2 y/R_(p) X-Tilt 2 2 1 +1 2r cos(θ) = 2/R_(p) 0 2 x/R_(p) Y- 3 5 2 −2 {square root over (6)}r² sin(2θ) = 2{square root over (6)}y/R_(p) ² 2{square root over (6)}x/R_(p) ² Astigmatism 2{square root over (6)} xy/R_(p) ² Defocus 4 4 0   0 {square root over (3)}(2r² − 1) = 4{square root over (3)}x/R_(p) ² 4{square root over (3)}y/R_(p) ² {square root over (3)}(2(x² + y²)/R_(p) ² − 1) X- 5 6 6 +2 {square root over (6)}r² cos(2θ) = 2{square root over (6)}x/R_(p) ² −2{square root over (6)}y/R_(p) ² Astigmatism {square root over (6)} (x² − y²)/R_(p) ² Y- 6 9 3 −3 {square root over (8)}r³ sin(3θ) = 6{square root over (8)} xy = 6{square root over (2)}(x² − y²)x = Trefoil {square root over (8)}(−y³ + 3x²y)/R_(p) ³ 6{square root over (8)} xy/R_(p) ³ {square root over (8)}(3x² − 3y) Y-Coma 7 7 3 −1 {square root over (8)}(3r³ − 2r) sin(θ) ={square root over (8)}(3x²y/R_(p) ³ + 3y³/R_(p) ³ − 2y/R_(p)) 6{square root over (8)}xy/R_(p) ³ $\sqrt{8}\left( {{3\frac{x^{2}}{R_{p}^{3}}} + {9\frac{y^{2}}{R_{p}^{3}}} - \frac{2}{R_{p}}} \right)$ X-Coma 8 8 3 +1 {square root over (8)}(3r³ − 2r) cos(θ) ={square root over (8)}(3x³/R_(p) ³ + 3 y²x/R_(p) ³ − 2x/R_(p)) $\sqrt{8}\left( {{3\frac{y^{2}}{R_{p}^{3}}} + {9\frac{x^{2}}{R_{p}^{3}}} - \frac{2}{R_{p}}} \right)$ $6\sqrt{8}\frac{xy}{R_{p}^{3}}$ X- Trefoil 9 10 3 +3 {square root over (8)}r³ cos(3θ) ={square root over (8)}(x³/R_(p) ³ − 3 y²x/R_(p) ³) $\sqrt{8}\frac{\left( {{3x^{2}} - {3y^{2}}} \right)}{R_{p}^{3}}$ ${- 6}\sqrt{8}\frac{yx}{R_{p}^{3}}$ Y- Quadrafoil 10 15 4 −4 {square root over (10)}r⁴ sin(4θ) = 4{square root over (10)} yx(x² − y²)/R_(p) ⁴ $4\sqrt{10}\frac{\left( {{3x^{2}y} - y^{3}} \right)}{R_{p}^{4}}$ $4\sqrt{10}\frac{\left( {x^{3} - {3xy^{2}}} \right)}{R_{p}^{4}}$ Y- 11 13 4 −2 {square root over (10)}(4r⁴ − 2{square root over (10)}(4 y³/R_(p) ⁴ + 2{square root over (10)}(4x³/R_(p) ⁴ + Secondary 3r²) sin(2θ) = 12 yx²/R_(p) ⁴ − 12 yx²/R_(p) ⁴ − Astigmatism {square root over (10)}(8 x³y/R_(p) ⁴ + 3 y/R_(p) ²) 3 x/R_(p) ²) 8 y³x/R_(p) ⁴ − 6 xy/R_(p) ²) Primary 12 11 4   0 {square root over (5)}(6r⁴ − 6r² + 1) = {square root over (5)}(24 x³/R_(p) ⁴ + {square root over (5)}(24 y³/R_(p) ⁴ + Spherical {square root over (5)}(6(x⁴ + y⁴ + 2x² − 24 y²x/R_(p) ⁴ − 24 x²y/R_(p) ⁴ − 6 x²/R_(p) ² − 6y²/R_(p) ² + 1) 12 x/R_(p) ²) 12 y/R_(p) ²) X- Secondary 13 12 4 +2 {square root over (10)}(4r⁴ − 3r²) cos(2θ) = $\sqrt{10}\left( {{16\frac{x^{3}}{R_{p}^{4}}} - {6\frac{x}{R_{p}^{2}}}} \right)$ $\sqrt{10}\left( {{6\frac{y}{R_{p}^{2}}} - {16\frac{y^{3}}{R_{p}^{4}}}} \right)$ Astigmatism {square root over (10)}[(4x⁴ − 4y⁴)/R_(p) ⁴ − (3x² − 3y²)/R_(p) ²] X- Quadrafoil 14 14 4 +4 {square root over (10)}r⁴ cos(4θ) ={square root over (10)}(x⁴ + y⁴ − 6y²x²)/R_(p) ⁴ $\sqrt{10}\left( {{4\frac{x^{3}}{R_{p}^{4}}} - {12\frac{{xy}^{2}}{R_{p}^{4}}}} \right)$ $\sqrt{10}\left( {{4\frac{y^{3}}{R_{p}^{4}}} - {12\frac{{yx}^{2}}{R_{p}^{4}}}} \right)$

In addition, Zernike coefficients may be derived in different ways, as is well known in the art. For example, an astigmatism prescription may be converted to a Zernike phase shift function and vice-versa. The corresponding coefficients of concern here are: S: the spherical dioptric power 1026, C: cylindrical dioptric power 1030, ϑ the cylinder axis angle 1032, and R_(P): the pupil radius. In this example, the relation to the relevant Zernike coefficients 1036 may be given through the following relations:

Plus-cylindrical notation:

$\begin{matrix} {{S = {{- \frac{4\sqrt{3}c_{2}^{0}}{R_{P}^{2}}} - \frac{2\sqrt{6}\sqrt{\left( c_{2}^{- 2} \right)^{2} + \left( c_{2}^{2} \right)^{2}}}{R_{P}^{2}}}},} \\ {{C = {- \frac{4\sqrt{6}\sqrt{\left( c_{2}^{- 2} \right)^{2} + \left( c_{2}^{2} \right)^{2}}}{R_{P}^{2}}}},} \\ {\vartheta = {\frac{1}{2}{\tan^{- 1}\left( \frac{c_{2}^{- 2}}{c_{2}^{2}} \right)}}} \end{matrix}$ Minus-cylindrical notation:

$\begin{matrix} {{S = {{- \frac{4\sqrt{3}c_{2}^{0}}{R_{P}^{2}}} + \frac{2\sqrt{6}\sqrt{\left( c_{2}^{- 2} \right)^{2} + \left( c_{2}^{2} \right)^{2}}}{R_{P}^{2}}}},} \\ {{C = {- \frac{4\sqrt{6}\sqrt{\left( c_{2}^{- 2} \right)^{2} + \left( c_{2}^{2} \right)^{2}}}{R_{P}^{2}}}},} \\ {\vartheta = {{\frac{1}{2}\tan^{- 1}\left( \frac{c_{2}^{- 2}}{c_{2}^{2}} \right)} + \frac{\pi}{2}}} \end{matrix}$ Zernike coefficients for plus and minus cases:

$\begin{matrix} {c_{2}^{- 2} = \frac{R_{P}^{2}C{\sin\left( {2\vartheta} \right)}}{4\sqrt{6}}} \\ {c_{2}^{0} = \frac{- {R_{P}^{2}\left( {S + {C/2}} \right)}}{4\sqrt{3}}} \\ {c_{2}^{2} = {\frac{R_{P}^{2}C{\cos\left( {2\vartheta} \right)}}{4\sqrt{6}}.}} \end{matrix}$

Thus, for example, considering an exemplary astigmatism prescription (no spherical dioptric power, cylindrical dioptric power of 2 and cylinder axis of 30°), an exemplary pupil radius of 2.5 mm and a refractive index of the phase element is taken to be equal to 1.55, an exemplary height profile function h(x, y) will be derived below.

In this example, the corresponding Zernike coefficients are given by:

$\begin{matrix} {c_{2}^{- 2} = {\frac{R_{P}^{2}C{\sin\left( {2\vartheta} \right)}}{4\sqrt{6}} = {0.17678R_{P}^{2}}}} \\ {c_{2}^{0} = {\frac{- {R_{P}^{2}\left( {S + {C/2}} \right)}}{4\sqrt{3}} = {- 0.14434R_{P}^{2}}}} \\ {c_{2}^{2} = {\frac{R_{P}^{2}C{\cos\left( {2\vartheta} \right)}}{4\sqrt{6}} = {0.10206{R_{P}^{2}.}}}} \end{matrix}$

The corresponding normal vector function for the second surface is given by:

$\begin{matrix} {{{\hat{n}\left( {x,y} \right)} = \frac{\left\langle {{- {\sum{C_{i}\frac{\partial{Z_{i}\left( {x,y} \right)}}{\partial x}}}},{- {\sum{C_{i}\frac{\partial{Z_{i}\left( {x,y} \right)}}{\partial y}}}},\left( {n - 1} \right)} \right\rangle}{\sqrt{\left( {\sum{C_{i}\frac{\partial{Z_{i}\left( {x,y} \right)}}{\partial x}}} \right)^{2} + \left( {\sum{C_{i}\frac{\partial{z_{i}\left( {x,y} \right)}}{\partial y}}} \right)^{2} + \left( {n - 1} \right)^{2}}}};} \\ {{\hat{n}\left( {x,y} \right)} = {\frac{\left\langle {{- \left( {{0.86604y} - {0.50003x}} \right)},{- \left( {{0.86604x} - {1.5y}} \right)},0.55} \right\rangle}{\sqrt{\left( {{0.86604y} - {0.50003x}} \right)^{2} + \left( {{0.86604x} - {1.5y}} \right)^{2} + (0.55)^{2}}}.}} \end{matrix}$

The corresponding height or thickness function h(x, y) is given by:

$\begin{matrix} {{{\varphi\left( {x,y} \right)} = {{{\mathbb{Z}}\left( {x,y} \right)} = {\frac{2\pi}{\lambda_{0}}{\sum{C_{i}{Z_{i}\left( {x,y} \right)}}}}}};} \\ {{= {\frac{2\pi}{\lambda_{0}}\left( {{0.86604{xy}} - {0.50001\left( {x^{2} + y^{2}} \right)} + {0.25R_{P}^{2}} + {0.24999\left( {x^{2} - y^{2}} \right)}} \right)}};} \\ {{= {\frac{2\pi}{\lambda_{0}}\left( {{0.86604{xy}} - {0.25002x^{2}} - {0.75y^{2}} + {0.25R_{P}^{2}}} \right)}};} \end{matrix}$ and thus:

${h\left( {x,y} \right)} = {{\frac{\lambda_{0}}{2\pi}\frac{\varphi\left( {x,y} \right)}{n - n_{0}}} = {\left( {{1.5746{xy}} - {0.45458x^{2}} - {1.3636y^{2}} + {0.25R_{P}^{2}}} \right).}}$

As an example, FIGS. 33A and 33B show plots of the resulting height function h(x, y) computed above (in meters).

Thus, going back to FIG. 31B, to build or construct functions h(x, y) and the associated normal vector functions, starting at sub-step 3106, for each coefficient in Zernike coefficient(s) 1036, at sub-step 3108 the corresponding contribution to the height or thickness function (and normal vector functions) for this Zernike term (coefficient and function) is added, as was done in the example above. This is done until all non-zero Zernike coefficients 1036 have been iterated upon at step 3110.

Going back to FIG. 31A, once the geometry or shape of a phase element 2500 is known, at sub-step 3104, it is “virtually” positioned at the correct location and with the correct orientation. In the examples given herein, first surface 2700 of phase element 2500 is always considered to be parallel to the (x, y) plane, but more general orientations may be considered. Computations with non-parallel surfaces may take into account the required rotational and translational transformations as necessary.

In addition, as mentioned above, phase element may be located at a given distance from the user's eye (so to model a real external lens) or it may be located (via the position of the first surface 2700) at pupil plane 1406 (to model an eye aberration).

Going back to FIG. 30 , the rest of sub-step 3002 continues as discussed above in the context of FIG. 13 (i.e. sub-steps 1301 to 1312), again depending on whether ray-tracing is done using the virtual image plane or the retinal plane.

Going back to FIG. 29 , the method proceeds with steps 1104 and 1106 as explained above for method 1100. However, at this point, ray 1414 (going through center location 1410 of the corresponding optical unit and originating from pixel location 1402) will encounter and be refracted by a phase element 2500, either before or after encountering pupil plane 1406.

Thus, step 2908, as illustrated in FIG. 35 , is where rays are propagated and interacted with each phase element 2500. Notably, in this embodiment, more than one phase element 2500 may be used simultaneously, as illustrated schematically in FIG. 44 . There, as an example, we see three different phase element 2500 located and oriented differently between the pixel display 106 and LFSL 108 on one side, and the pupil plane 1406 on the other side. These phase elements 2500 may or may not have the same optical properties, as required. In addition, this exemplary embodiment of step 2908 further includes the possibility of having a phase element 2500 at the pupil plane itself, so as to model an aberrated eye, or direct compensation therefor.

Going back to FIG. 35 , at sub-step 3502, starting from optical element center location 1410, ray 1414 is propagated in a straight line towards the eye until it either encounters an element 2500 in its path (via first surface 2700), or until it reaches the pupil plane 1406.

At sub-step 3503, a check is made to see if the ray as reached the pupil plane. If not, then it has encountered a first surface 2700 of a phase element 2500 in its path. At sub-step 3504 the refracted ray obtained from this phase element 2500 is computed.

This is schematically illustrated in FIG. 34 , where phase element 2500 is shown to be located somewhere between the eye and the light field display. Thus, in this example, first surface 2700 is shown being located at a distance 3450 from pupil plane 1406. In the case where the eye itself is being modelled, then distance 3450 would be zero and first surface 2700 would be at pupil plane 1406. In addition, ray 1414 is shown encountering a first surface 2700 of a phase element 2500 at input location 3404, which produces a corresponding refracted ray 3402 exiting from second surface 2702 at output location 3406. Thus, at sub-step 3504, both the origin and orientation of refracted ray 3402 exiting phase element 2500 is computed. It is this ray that is propagated.

The skilled technician will understand that, at sub-step 3504, different methods or algorithms may be used to compute the origin and orientation of refracted ray 3402 coming out of the second surface 2702 of phase element 2500. In general, this involves computing a refraction event via electromagnetic wave propagation equations or Snell's law for each surface. However, in some embodiments, different approximations may be used to minimize the computational load of such computations.

As shown schematically in FIG. 36A, in the currently discussed embodiment (where first surface 2700 is flat and the second surface 2702 has a changing profile h(x, y)), incoming ray 1414 will intersect the first surface 2700 at input location 3404 (herein corresponding to x, y coordinates (x₀, y₀), be refracted as ray 3602, which will travel within phase element 2500 until it reaches second surface 2702 at output location 3406 (herein corresponding to x, y coordinates (x₁, y₁)), where it will once again be refracted. The location in z of output location 3406 is given by the height or thickness profile function at this point, i.e. h(x₁, y₁). Refracted ray 3402 is thus the refracted ray or ray vector exiting phase element 2500 at output point 3406.

In this exemplary embodiment, the normal vector {circumflex over (n)} to first surface 2700 is by definition parallel to the z-axis, while the normal vectors of the second surface 2702 are defined as a function of x, y coordinates: {circumflex over (n)}(x, y), as was discussed above. Thus, in FIG. 36A, the normal vector at second surface 2702 is given by {circumflex over (n)}(x₁, y₁).

As mentioned above, different methods or techniques may be used to efficiently compute the refraction events illustrated in FIG. 36A. Some examples are discussed below.

In the most computationally demanding example, ray 3602 is projected forward within phase element 2500 using small propagation steps until the height or thickness traversed is equal or larger than the height or thickness at the current iteration, at which point the second refraction event is computed. However, this requires the use of small steps and thus many iterations may be required.

In one exemplary approximation, illustrated schematically in FIG. 36B, for a substantially thin phase element, output location 3406 may be approximated (e.g. thin lens approximation) as being the same as input location 3404 (i.e. (x₀, y₀). Thus, the corresponding z-coordinate and normal vectors are given by h(x₀, y₀) and normal vector {circumflex over (n)}(x₀, y₀).

In another example, as illustrated in FIG. 36C, ray 3602 may be propagated inside phase element 2500 a distance equal to the height or thickness at the first intersection position h(x₀, y₀) and then the normal vector at the new coordinates (x₃, y₃) may be used to compute the second refraction event.

In addition, as mentioned above, in some embodiments, the values of functions h(x, y) and {circumflex over (n)}(x, y) may be pre-computed and stored in an array-like data structure, for example in a texture file. An example is given schematically in FIGS. 37A to 37C. Thus, the exemplary phase element 2500 of FIG. 37A (herein shown to have a rectangular shape as an example only to better illustrate) may be stored as a pre-computed array of height or thickness values herein represented as digitized phase element 4402 in FIG. 37B. Thus, each portion (i, j) (i.e. array index number) corresponds to a given height value h(i, j) and normal vector {circumflex over (n)}(i, j). The resolution of phase element 4402 may be changed as required by computational considerations, thus varying the corresponding size of each portion (dx, dy). FIG. 37C shows a cross-sectional view of FIG. 37B superimposed on the surfaces 2700 and 2702 of the exemplary phase element of FIG. 37A. In the case where the values are stored within a texture data structure, each portion corresponds to a texel (texture pixel) or similar and thus a corresponding (x, y) coordinate (for example input location (x, y) located within a given portion may be converted to the appropriate index (i, j) to access the corresponding height and normal vector).

In some embodiments, interpolated values for the height and normal vector may be generated as a function of the relative position of input location (x, y) with respect to the neighboring texel locations. This may include for example linear interpolation but other interpolation methods may also be used.

In some embodiments, the pre-computed set of values h(i, j) may be sorted from smallest to largest, and then checked iteratively while propagating ray 3602 until the height or thickness at the newly obtained (x, y) coordinates is within a given tolerance interval.

Notably, as illustrated in FIG. 34 , refracted ray 3402 will be perceived by the user, via projection 3408, as originating from a “virtual” optical element center 3410 (e.g. “virtual” lenslet or microlens center position) and corresponding “virtual” pixel 3412. As will be discussed further below, in some embodiments, each set of “virtual” pixel 3412 and associated “virtual” optical element center 3410 computed from each pixel in the pixel display, as discussed above, may be computed once and kept or stored in memory for later iterations.

Going back to FIG. 35 , once a first refracted ray 3402 has been computed, the ray is propagated once more via sub-step 3502 until it encounters another phase element 2500 or until it reaches the pupil. Thus, in the case where multiple phase elements 2500 are present between the pupil plane 1406 and the LFSL 108, sub-steps 3502, 3503 and 3504 may be repeated as many times as necessary. Each time, the last refracted ray 3402 is used to compute the next refraction event. It is to be understood that the refracted ray 3402 that reaches the pupil plane 1406 is the one resulting from the last refraction event (e.g. last execution of sub-step 3504).

When, at sub-step 3503, refracted ray 3402 (or ray 1414 if no phase element 2500 has been encountered thus far) reaches the pupil plane 1406, the method proceeds to sub-step 3506 where a check is made that the intersection point of ray (1414 or 3502) on pupil plane 1406 is still within the pupil entrance, just like in step 1108 of method 1100 for example. If this is not the case, then the user cannot see the ray and step 2908 exits with a NO condition, so that method 2900 proceeds to step 1110 in FIG. 29 . In the opposite case, then a final check is made at sub-step 3508 to see if a phase element 2500 has been placed at the pupil plane itself (to model an aberrated or visually compensated eye). If so, then sub-step 3510 proceeds just like sub-step 3504 and the new refracted ray 3402 is computed, just as discussed above. If no phase element was present at the pupil plane at sub-step 3508, or when the new refracted ray 3402 is computed at sub-step 3510, step 2908 exits with the YES condition, so that in method 2900 of FIG. 29 , the method then proceeds to step 2912.

There, the last computed refracted ray 3402 (and its projection 3408 and “virtual” optical element center location 3410) may be used to replace ray 1414 in step 2912, which is illustrated in FIG. 38 and comprises the same sub-steps as step 1112 of FIG. 16 , but with some changes, as will be discussed below.

Thus, step 2912 proceeds via sub-step 1601 as discussed above where the path changes depending on whether the ray-tracing is done on the virtual image plane or on the retinal plane. If it's the former, then sub-step 1602 is executed to identify the image portion on virtual image plane 1502, just like as described above in the context of FIG. 16 , but here it is the intersection point of projection vector 3408 on virtual image plane 1502 (instead of vector 1504) which provides the corresponding image portion 1506 in FIG. 15 .

In contrast, if at sub-step 1601 ray-tracing is done on the retinal plane, then sub-step 3808 is executed. Here, instead of using the real optical element center location 1410 to draw vector 1718 through pupil center 1022, it is instead the “virtual” optical element center location 3410 that is used instead. Thus, the ray-tracing steps illustrated in FIGS. 17A to 17D apply equally here as well, with the exception that point 1410 in FIG. 17A is replaced by point 3410. Thus, in FIG. 38 , sub-steps 1610 to 1614 proceed as explained above to identify the image portion 1724.

Going back to FIG. 29 , once step 2912 is finished and the corresponding image portion identified, steps 1114, 1116 and 1118 proceed just like in method 1100.

With reference to FIGS. 39 to 43 , and in accordance with one exemplary embodiment, a modified version of method 2900 for correcting or compensating for HOAs but also operable to render multiple images or image portions on multiple adjusted distinct image planes simultaneously, generally referred to using the numeral 3900, will now be discussed.

Method 3900 is to method 2900 what the steps discussed above and illustrated in FIGS. 20A to 23B were for method 1100: it allows the light field display to render multiple optotypes or images simultaneously at different correction levels. In this exemplary embodiment, method 3900 as discussed below further includes the possibility of using multiple phase elements, but here wherein one distinct phase element is used for each image or optotype being rendered, and wherein each element has its own optical properties and thus simulates a different set of HOAs (or aberrations in general). Each set of HOAs (or aberrations in general) may also be defined via its own set of Zernike coefficients 1036. Thus, multiple corrections to HOAs (or aberrations in general) may be rendered simultaneously in different image portions or distinct optotypes side-by side. Usually, each set of HOAs may also be used in conjunction with a corresponding spherical dioptric power 1026, which will be used to define either a corresponding virtual image plane 1502 or an eye focal plane 1716, depending on the ray-tracing path chosen.

Thus, in FIG. 39 , method 3900 begins with step 3902, which mirrors in part step 2002 of FIG. 20A. However, as illustrated in FIG. 40A, step 3902 further includes, prior to sub-step 2400, sub-step 4002 where multiple phase elements may be constructed (for example one for each optotype being rendering simultaneously). Sub-step 4002 is illustrated in the process flow diagram of FIG. 41 , where at sub-step 4102 a corresponding correction for a given set of HOAs (or aberrations in general) (for a given image portion or optotype) is selected (and its corresponding set of Zernike coefficients 1036). Then this set of Zernike coefficients is used to generate and place a corresponding phase element 2500 in sub-step 3102 and 3104 as discussed above. Then at sub-step 4104, if other sets of HOAs (or aberrations in general) are to be included, the method goes back iteratively to sub-step 4002 as many times as necessary, until all corresponding phase elements have been generated. Going back to FIG. 40A, after sub-step 4002, steps 2004, 2102 and 2300 proceed just as described above for each virtual image plane (sub-step 2102) or eye focal plane (sub-steps 2300).

Going back to FIG. 39 , the iteration over each pixel proceed as before, via steps 1104 and 1106 discussed above. Then at step 3908, which is further illustrated in FIG. 40B, if at sub-step 2004 ray-tracing is done using virtual image planes, then the method proceeds to sub-step 4004, which is illustrated in FIG. 42 , while if ray-tracing proceeds on the retinal plane, the method continues proceeds with sub-step 4006, which is illustrated in FIG. 43 .

In the former case (virtual image planes), in some embodiments, sub-step 3810 is similar to step 2104 of FIG. 21B, but includes at sub-step 4210 a loop not only over corresponding spherical dioptric power 1026 values but also over associated (if any) HOAs and their corresponding phase element. So, once a combination of virtual image plane and phase element is selected, the method uses these to execute sub-step 3208, as described above. If, as computed within sub-step 3208, the ray (either ray 1414 or refracted ray 3402) does not reach inside the user pupil, then the method goes to sub-step 1110, described above. In contrast, if the raydoes reach within the pupil, then the method proceeds with sub-step 1602 as discussed above, wherein the image portion on the virtual image plane is identified. Then the color channel of this image portion is identified at sub-step 2112, as before. Then, at sub-step 2114, as discussed before, if the channel is transparent, then a check is made if all combinations of virtual image planes and phase elements have been used at sub-step 4216, if not then the method goes back to sub-step 4210. If, at sub-step 4216 all combinations have been used and no non-transparent color channel has been encountered, then the method goes to sub-step 2118, as discussed above, so that the color is set to black (or another background color).

Going back to FIG. 40 , if at sub-step 2004 ray-tracing is selected to be done on retina plane 1702 instead, then sub-step 4006 is executed. Flow process diagram of FIG. 43 illustrated one embodiment of sub-step 4006. Similarly, to sub-step 4004, a loop is executed over all sets of spherical dioptric power 1026 values but also over corresponding associated (if any) HOAs and their corresponding phase element at sub-step 4310. Then, sub-step 3208 is executed just as described above, and thus includes refraction by the corresponding phase element 2500 but also a check on whether the ray (either ray 1414 or refracted ray 3402, depending on the location of this particular phase element 2500) hits inside the pupil opening. If not, then step 1110 is executed as discussed above. If it does hit within the pupil opening, then steps 3808, 1610, 1612 and 1614 are executed one after the other as discussed above in the context of FIG. 38 , where the “virtual” optical unit center location 3410 is used instead of the optical unit center location 1410 to draw vector 1718. Once the correct image portion has been identified at sub-step 1614, then step 2308 is executed to assign the corresponding colour channel to pixel 1402. Then sub-step 2310 is executed where, as explained above, the colour channel is checked to see if it has been defined as transparent. If it is not, then the colour value is retained and sub-step 4006 is exited. If it is transparent, then the step 4312 is executed to see if all of the corresponding sets of eye focal planes and phase elements have been iterated upon. If not, then the method goes back to sub-step 4310 with the next eye focal plane and phase element. If no phase element and eye focal plane is left, then sub-step 2825 is executed as described before, so that the corresponding colour channel is set to black (or another background colour), before exiting sub-step 4006.

Going back to FIG. 39 , once step 3908 is finished, at step 1114 (described above) colour channel is assigned to pixel 1402. Steps 1116 (check if all pixels or sub-pixels have been processed) and 1118 (wherein the image is rendered via pixel display 110) proceed as explained before.

With reference to FIG. 45 , and in accordance with one exemplary embodiment, a unified light field rendering method incorporating phase elements, generally referred to using the numeral 4500, will now be described. Generally, this method is used to project a scene from a light field display using the “virtual” pixel and optical unit center (i.e. like “virtual” pixel 3412 and an associated “virtual” optical unit center 3410 discussed above and with reference to FIG. 34 ).

The method is to be implemented without the need of transforming the projected scene to an image plane or projection on retina. The role of the light-field algorithm, here, is to sample the scene/virtual object using the generated rays after considering eye aberrations using the phase element. The retinal spot size and spacing control is to be handled with optics design to control the generated beam characteristics.

Thus, at step 4502, a virtual object/scene is set at a position regardless of the refractive error (e.g. aberration) of the user's eye. Then at step 4504, an aberration-free eye is assumed to perceive the reconstructed scene without defects. At step 4506, any refractive error that aberrates the eye and needs to be corrected via light field rendering is lumped into one or more phase elements 2500, as discussed above, each element having the required position and orientation.

At step 4508, for each pixel 1402 in pixel display 108, rays are propagated through the corresponding optical element center 1410 and refraction events from the one or more phase elements 2500 are computed to derive a final refracted ray (i.e. ray 3402 in FIG. 34 ).

At step 4510, each ray 3402 is shifted or translated to the nodal point of the eye lens (e.g. to cross pupil center location 1022 for example), as illustrated schematically in FIG. 46 . There, we see a ray 3402 (with a corresponding virtual optical unit center location 3410 and associated virtual pixel location 3412) being shifted to intersect pupil center location 1022, thus giving ray 4602. By tracing back, another set of virtual pixel location 4604 and virtual optical unit center location 4606 may be found. In FIG. 46 , D_(LE) is the distance between the LFSL 106 and pupil plane 1406, D is the distance between the pupil and the retina (i.e. eye depth value 1034) which is assumed to be equal to the focal length (FL) of the aberration-free eye, and D_(PE) is the distance between the pupil plane 1406 and the pixel display 106.

At step 4512, raytracing may be done, using different methods, such as method 1100 for example, by using virtual pixel location 4604 and virtual optical unit center location 4606 in place of the real pixel location and optical unit center location (thus using ray 4602 instead of ray 1414). Thus, ray 4602 may, for example, by projected backwards to a scene or object 4608 (i.e. for example like ray 1504 on virtual image plane 1502) to identify an corresponding image portion to be assigned to the corresponding pixel.

Below are some examples of how method 4500 may be applied to different situations:

Multiple Elements Treatment

For multiple elements with spatial offsets, method 4500 may be applied in a recursive manner or by full tracing of the beams generated by the lenslets-pixels combinations starting from the element closer to the lenslet (D_(PO1) to D_(PO3) in FIG. 44 ) and ending by the eye pupil, taking into account the required space transformation. For recursive implementation, the lenslets and pixel coordinates are updated after each iteration then the refracted phase element is removed for the following iteration. At the end of tracing, the virtual positions of lenslets and pixels are calculated, and light field raymarching is applied. For each phase element done with space transformation, inverse transformation can be done or the whole system is transformed including the virtual object all the way through.

Multi-Prescription Treatment

For a multi-prescription treatment, the scene objects/charts are sorted in space in the positions and scales they are desired to be viewed. For each prescription, a phase element is constructed, and virtual positions are calculated for the lenslets and pixels. Ray-tracing is done for each prescription/scene object combination. The pixels of rays that hit the prescription region being calculated are excluded from the following prescription calculation.

Progressive Depth Treatment

Elements like progressive lenses have different focal lengths at different regions on the lens element. For this kind of application, the phase element is constructed to correspond to refraction/phase profile of these elements then the regular procedure is applied as explained in steps 4502 to 4512 above.

While the present disclosure describes various embodiments for illustrative purposes, such description is not intended to be limited to such embodiments. On the contrary, the applicant's teachings described and illustrated herein encompass various alternatives, modifications, and equivalents, without departing from the embodiments, the general scope of which is defined in the appended claims. Except to the extent necessary or inherent in the processes themselves, no particular order to steps or stages of methods or processes described in this disclosure is intended or implied. In many cases the order of process steps may be varied without changing the purpose, effect, or import of the methods described.

Information as herein shown and described in detail is fully capable of attaining the above-described object of the present disclosure, the presently preferred embodiment of the present disclosure, and is, thus, representative of the subject matter which is broadly contemplated by the present disclosure. The scope of the present disclosure fully encompasses other embodiments which may become apparent to those skilled in the art, and is to be limited, accordingly, by nothing other than the appended claims, wherein any reference to an element being made in the singular is not intended to mean “one and only one” unless explicitly so stated, but rather “one or more.” All structural and functional equivalents to the elements of the above-described preferred embodiment and additional embodiments as regarded by those of ordinary skill in the art are hereby expressly incorporated by reference and are intended to be encompassed by the present claims. Moreover, no requirement exists for a system or method to address each and every problem sought to be resolved by the present disclosure, for such to be encompassed by the present claims. Furthermore, no element, component, or method step in the present disclosure is intended to be dedicated to the public regardless of whether the element, component, or method step is explicitly recited in the claims. However, that various changes and modifications in form, material, work-piece, and fabrication material detail may be made, without departing from the spirit and scope of the present disclosure, as set forth in the appended claims, as may be apparent to those of ordinary skill in the art, are also encompassed by the disclosure. 

What is claimed is:
 1. A computer-implemented method, automatically implemented by one or more digital processors, to automatically adjust user perception of an input image to be rendered via a set of pixels and a corresponding array of light field shaping elements (LFSE), wherein the user perception is affected by a higher order aberration, the method comprising: receiving as input one or more higher order aberration parameters digitally defining the higher order aberration; for each given pixel of at least some of said pixels, digitally identifying an adjusted image plane location corresponding thereto given a corresponding LFSE corresponding thereto and given said one or more higher order aberration parameters; digitally associating with said adjusted image plane location an adjusted image value designated for said adjusted image location; and rendering for each said given pixel said adjusted image value associated therewith.
 2. The computer-implemented method of claim 1, further comprising for said each given pixel: digitally projecting an adjusted image ray trace from said given pixel towards a user pupil location, given said corresponding LFSE corresponding thereto, to intersect a virtual phase element digitally defined between said given pixel and the user's retina to virtually produce a designated optical phase shift corresponding at least in part to the higher order aberration; digitally computing a corresponding refracted ray trace resulting from intersection of said adjusted image ray trace with said virtual phase element; and digitally projecting said refracted ray trace to intersect said adjusted image plane location on said designated adjusted image plane.
 3. The computer-implemented method of claim 2, wherein said virtual phase element is defined by a variable dioptric power profile, and wherein said refracted ray trace is refracted as a function of a given dioptric power defined for an intersection portion of said virtual phase element intersected by said adjusted image ray trace.
 4. The computer-implemented method of claim 3, wherein said variable dioptric power is profiled to simulate a progressive or bifocal lens.
 5. The computer-implemented method of claim 3, wherein said virtual phase element is defined by a height function characterizing said variable dioptric power profile.
 6. The computer-implemented method of claim 5, wherein said height function is derived from a set of Zernike coefficients preset so as to reproduce said designated optical phase shift.
 7. The computer-implemented method of claim 5, wherein said virtual phase element is further defined by a normal vector function derived from said height function that characterizes an orientation of an output surface of said virtual phase element with respect to any location on an input surface of said virtual phase element; wherein said height function and said normal vector function are represented in memory in a texture data object comprising, respectively, a set of height values and associated normal vector values of said output surface corresponding to a plurality of surface locations on said input surface; wherein using said height function and said normal vector function involves fetching a height value and associated normal vector value corresponding to the surface location of said plurality of surface locations that is closest to said input location.
 8. The computer-implemented method of claim 2, wherein said virtual phase element is defined by an input surface and an output surface, and wherein said refracted ray trace is digitally refracted to output said refracted ray trace by sequentially computing a first refraction event at said input surface and a second refraction event at said output surface, wherein at least one of said input surface or said output surface has a variable profile.
 9. The computer-implemented method of claim 2, wherein said input image comprises distinct image portions, each image portion being associated with a distinct virtual phase element, and wherein said adjusted image ray trace for said distinct image portions intersects a corresponding one of said distinct virtual phase element to be processed accordingly.
 10. The computer-implemented method of claim 2, wherein said virtual phase element comprises a plurality of virtual phase elements, and wherein said adjusted image ray trace is sequentially processed to intersect sequential ones of said virtual phase elements in its path.
 11. The computer-implemented method of claim 2, wherein said virtual phase element is virtually located at a location between said given pixel and a user pupil location so as to at least partially simulate an intervening corrective lens defined to produce said designated phase shift in at least partially addressing the optical aberration.
 12. The computer-implemented method of claim 2, wherein said virtual phase element is virtually located at a user's eye pupil or lens location so to simulate an optical correction to be experienced by the user after undergoing a corresponding eye surgery or intraocular implant.
 13. The computer-implemented method of claim 1, wherein said adjusted image plane comprises a pixel plane.
 14. The computer-implemented method of claim 1, further comprising receiving as input an eye focus parameter, wherein said adjusted image plane comprises a virtual image plane virtually positioned relative to the digital display at a designated minimum viewing distance designated such that said perceptively adjusted version of the input image is adjusted to accommodate said eye focus parameter.
 15. The computer-implemented method of claim 1, further comprising receiving as input an eye focus parameter, wherein said adjusted image plane is designated as a user retinal plane, and wherein said refracted ray trace is projected to intersect with said user retinal plane by further redirecting said refracted image ray trace given said pupil location in accordance with said eye focus parameter.
 16. A non-transitory computer-readable medium, comprising instructions for implementation by one or more digital processors to automatically adjust user perception of an input image to be rendered via a set of pixels and a corresponding array of light field shaping elements (LFSE), wherein the user perception is affected by a higher order aberration, by: for each given pixel of at least some of said pixels, digitally identifying an adjusted image plane location corresponding thereto given a corresponding LFSE corresponding thereto and given one or more higher order aberration parameters defining the higher order aberration; digitally associating with said adjusted image plane location an adjusted image value designated for said adjusted image location; and rendering for each said given pixel said adjusted image value associated therewith.
 17. The non-transitory computer-readable medium of claim 16, wherein the instructions further comprise, for said each given pixel: digitally projecting an adjusted image ray trace from said given pixel towards a user pupil location, given said corresponding LFSE corresponding thereto, to intersect a virtual phase element digitally defined between said given pixel and the user's retina to virtually produce a designated optical phase shift corresponding at least in part to the higher order aberration; digitally computing a corresponding refracted ray trace resulting from intersection of said adjusted image ray trace with said virtual phase element; and digitally projecting said refracted ray trace to intersect said adjusted image plane location on said designated adjusted image plane.
 18. The non-transitory computer-readable medium of claim 17, wherein said virtual phase element is defined by a variable dioptric power profile, and wherein said refracted ray trace is refracted as a function of a given dioptric power defined for an intersection portion of said virtual phase element intersected by said adjusted image ray trace.
 19. The non-transitory computer-readable medium of claim 17, wherein said virtual phase element is defined by an input surface and an output surface, and wherein said refracted ray trace is digitally refracted to output said refracted ray trace by sequentially computing a first refraction event at said input surface and a second refraction event at said output surface, wherein at least one of said input surface or said output surface has a variable profile.
 20. The non-transitory computer-readable medium of claim 17, wherein said input image comprises distinct image portions, each image portion being associated with a distinct virtual phase element, and wherein said adjusted image ray trace for said distinct image portions intersects a corresponding one of said distinct virtual phase element to be processed accordingly.
 21. The non-transitory computer-readable medium of claim 17, wherein said virtual phase element is virtually located at a location between said given pixel and a user pupil location so as to at least partially simulate an intervening corrective lens defined to produce said designated phase shift in at least partially addressing the optical aberration.
 22. The non-transitory computer-readable medium of claim 17, wherein said virtual phase element is virtually located at a user's eye pupil or lens location so to simulate an optical correction to be experienced by the user after undergoing a corresponding eye surgery or intraocular implant. 